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Abstract 


In this paper, we study the holographic Renyi entropy of a large interval on a circle 
at high temperature for the two-dimensional CFT dual to pure AdSa gravity. In the field 
theory, the Renyi entropy is encoded in the CFT partition function on n-sheeted torus 
connected with each other by a large branch cut. As proposed in [^, the effective way 
to read the entropy in the large interval limit is to insert a complete set of state bases 
of the twist sector at the branch cut. Then the calculation transforms into an expansion 
of four-point functions in the twist sector with respect to e ^. By using the operator 
product expansion of the twist operators at the branch points, we read the first few terms 
of the Renyi entropy, including the leading and next-to-leading contributions in the large 
central charge limit. Moreover, we show that the leading contribution is actually captured 
by the twist vacuum module. In this case by the Ward identity the four-point functions 
can be derived from the correlation function of four twist operators, which is related to 
double interval entanglement entropy. Holographically, we apply the recipe in 20 and 26 


to compute the classical Renyi entropy and its 1-loop quantum correction, after imposing 
a new set of monodromy conditions. The holographic classical result matches exactly with 
the leading contribution in the field theory up to and U, while the holographical 

1-loop contribution is in exact agreement with next-to-leading results in field theory up to 
and U as well. 
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1 Introduction 


The entanglement entropy is an important notion in a quantum many-body system [^[^. Not 
only could it be used to measure the effective degrees of freedom in the system, but it could 
also be taken as a quantum order parameter, among its various applications. It is defined 
as follows. Let A be a subsystem, and then the reduced density matrix of A is obtained by 
tracing out the degrees of freedom of its complements'^ 

PA = TrA=/9, (1.1) 

where p is the density matrix of the whole system. Then the entanglement entropy is defined 
to be the von Neumann entropy of the reduced density matrix 

Sa =-TrpAlogpA- (1-2) 

Furthermore for pure state p =\ ip){'4’ I; the entanglement entropy of the subsystem is equal 
to the one of its complementary part 

5a = 5ac, (1.3) 

but for a thermal state the equality breaks down, 

5a ^ 5a= (1.4) 

because of the thermal effect. It is convenient to calculate the entanglement entropy from the 
Renyi entropy, which is defined to be 

Sn = -^ logTrp)|. (1.5) 

The entanglement entropy can be read from 

See = lim Sn, (1.6) 

n—¥l 

if the limit n —)• 1 is well defined. 

In quantum field theory, the entanglement entropy and Renyi entropy are hard to compute 
because there are an infinite number of degrees of freedom. In this case, the entanglement 
entropy is defined with respect to a spatial submanifold at a fixed time. By using the replica 
trick the Renyi entropy can be transformed into the partition function of n copies of 
field theory with the fields being identified at the submanifold. It is usually a formidable 
task to compute this partition function for a general field theory. Even for two-dimensional 
(2D) conformal field theory (GET), which is expected to give more analytic results due to 
the existence of infinite dimensional symmetries, the exact results are limited. For 2D CFT, 
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the Renyi entropy is generally related to the partition function on a higher genus Riemann 
surface. Besides a few universal results determined by the conformal symmetries [^, only the 
partition functions of a free boson and fermion on a higher genus Riemann surface have been 
known U 0 d i 


However, it is possible to expand the partition function with respect to some modular 
parameters for a general CFT in some cases. The two simplest nontrivial examples are the 
case of two intervals on a complex plane and the case of one interval on a torus. Generally, to 
calculate the partition function, one can cut the Riemann surface at some cycles and insert a 
complete set of state bases such that the full Riemann surface changes into a surface without 
handle and hole and the computation transforms into a summation of multipoint correlation 
functions on a full complex plane. The key point is to find the nice way to cut open the Riemann 
surface such that the expansion series is well behaved. For a general genus-i^ Riemann surface. 


we can always choose g couples of A cycles and B cycles with a proper intersection 11 and cut 
the Riemann surface at certain cycles. Different choices on the cutting correspond to different 
ways of expanding the partition function. Even though by the modular invariance the different 
expansions should be equal to each other, their convergent rates are different. 

The simplest trivial example is the partition function on a torus. One may quantize the 
theory along the thermal direction or the spatial direction, which corresponds to inserting the 
complete bases along the spatial cycle or thermal cycle, and the partition function could be 
written as 


Z = 


E 




or Z = 


E 




(1.7) 


Because of the modular invariance, the two different calculations give the same answer. At 
a low temperature the quantization along the thermal direction leads to a better convergent 
series while at a high temperature the spatial quantization works better. 

For the double-interval case, if the intervals are short, one may take the operator product 
expansion(OPE) of the twist operators to compute the Renyi and entanglement entropy order 


by order with respect to a small cross ratio 12,13 . Actually, taking OPE is equivalent to 
inserting a complete set of bases at the cycles around the two intervals on every sheet. In 
the Riemann surface for the nth Renyi entropy of the double-interval case, there are n — 1 
independent couples of cycles denoted by Ai and Bj. As shown in Eig. [T]for the n = 3 case, 
there are two couples of independent cycles Ai, Bi,i = 1,2. In the small interval limit, we 
can take the OPE of the twist operators at the branch points of the hrst interval. This is 
equivalent to cutting and inserting a complete set of state bases at Ai cycles enclosing the 
interval. This expansion is well convergent for a small cross ratio. On the contrary, for a 
large cross ratio which means the intervals are large and the branch points of two intervals are 
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Figure 1: Canonical cycles for the double intervals 


close to each other, we must take the OPE of the second and third twist operators, which is 
equivalent to cutting the Riemann surface along the Bi cycles enclosing the branch points of 
separated intervals. 

For the case of a single interval on a torus, the Riemann surface is obtained by connecting 
n tori along the branch cut. When the interval is not very large, one may cut the spatial or 
thermal cycle and insert a complete set of bases to compute the partition function [14| , 
just as in the genus-1 case. Which cycles to cut depends on the temperature. But for a high 
temperature and a very large interval, the previous treatment is not good enough because 
the resulting expansion series is poorly behaved. Instead, it was proposed in that one 

should cut the cycle crossing the branch cut. This requires the insertion of a complete set of 
bases in the twist sector rather than the normal sector in the n-copied CFT. This proposal has 
been checked for the free compact and noncompact bosons, and it has been applied to prove 
the universal relation between the thermal entropy and the entanglement entropy. 

The AdS/CFT correspondence provides another way to compute the entanglement entropy 
in a CFT. For the Einstein gravity, it was first proposed by Ryu and Takayanagi |16||17| that 
the entanglement entropy could be holographically given by the area of a minimal surface in 
the bulk, which is homogeneous to A 


See 


Area(SA) 

4G 


( 1 . 8 ) 


The holographic entanglement entropy could be understood as a generalized gravitational en¬ 
tropy [18] |19j . In the higher dimension case, it is not clear if the holographic entanglement 
entropy gives precisely the entanglement entropy in the dual field theory. Nevertheless, for a 
2D CFT holographically dual to AdSs gravity, it has been proved that the holographic com¬ 


putation is correct in the semiclassical regime 20,21 . Therefore, the Renyi entropy provides 
a new window to study the AdS 3 /CFT 2 correspondence. 
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The AdS^/CFT 2 correspondence states that the quantum gravity in AdS'i spacetime is 
dual to a 2D CFT with a central charge 

31 


c = 


2G’ 


(1.9) 


and a sparse light spectrum [21| [23| , where G is the three-dimensional (3D) gravity coupling 
constant. Though a precise definition of AdSs quantum gravity, possibly a string theory, 
has not been well established, its semiclassical limit has been much studied. As the classical 
configurations in the AdSs gravity could be obtained as the quotients of the global AdSs by 
the subgroup of the isometry group SL{2,C), the path integral of semiclassical AdSs gravity 
could be defined in principle. On the other side, the explicit construction of dual CFT is 
not known. Nevertheless, the large central charge limit of the CFT, corresponding to the 
semiclassical gravity, is much simplified. Under this limit, only the vacuum module dominates 
the contribution to the CFT partition function [^. As a result, the partition function is 
universal in the sense that it is very much restricted by the conformal symmetry, and is 
independent of the explicit construction of the CFT. In this work, we are interested in the large 
central charge limit of the Renyi entropy of 2D CFT. From the AdS 3 /CFT 2 correspondence, 
the partition function of the Riemann surface in the CFT should be given by the partition 
function of the gravitational configuration ending on the Riemann surface. In the large central 
charge limit, the Renyi entropy can be decomposed into the terms proportional to c, c°, ..., 

which should correspond respectively to the classical, quantum 1-loop, 2-loop, etc. , parts of 


the gravitational partition function 13 


In the field theory side, for a genus-fif Riemann surface, we need to choose g Ai cycles to 
insert complete bases at each cycle such that the expansion converges fast. Under the large 
c limit, the dominant contribution to the partition function comes from the light primary 
states and their descendants. The heavy states give only nonperturbative corrections of 
order 0{e~^). Furthermore, among the light spectrums only the vacuum module gives the 
linear c order result and the other modules give only higher order corrections with respect to 
. Moreover, it turns out that even for the next-to-leading correction with respect to -, 


21 


the first few terms in the expansion are captured only by the vacuum module 25 [14| . 
The vacuum Verma module consists of a primary identity operator and its descendants which 
could be constructed by the stress tensors T{z) and T{z). In this work, we assume that for 
the field theory that is dual to the pure gravity, we only need to insert the vacuum module at 
each cycle. 

Because of the replica symmetry, we always deal with a CFT on an n-sheeted surface, which 
can be regarded as one CFT with n copies of the original field with the fields in different replicas 
being identified along the branch cut. When we combine the n-sheeted surface’s field into one 
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CFT, we call it n-copied CFTj^with nc central charge and denote it as CFT„. The original 
CFT is denoted as CFTi. If we do not consider the monodromy condition of the fields around 
the branch point, the n copies of the fields are decoupled so that we have just a tensor product 
of n copies of the fields. In this case, we call it the normal sector of CFT„. In contrast, if we 
consider the twist monodromy condition of the fields around the branch point, we get the twist 
sector of CFT,i. In both cases, we can classify the states by the irreducible representations 

n 

of its Virasoro algebra Vir^*), defined by the stress tensor T{z) = E T^^\z). Under Vir^*^ 

i=i 

the twist (normal) sector states can be decomposed into more than one irreducible module, 
and the one with the lowest conformal dimension is called the twist (trivial) vacuum module. 
Note that the twist (trivial) vacuum module has a different meaning from the vacuum module 
in the original CFT. For the partition function of a CFT on an re-sheeted Riemann surface 
resulting from the replica trick, we have two pictures to compute it. One is to regard it as the 
one of CFTi on the re-sheeted surface, and the other one is to regard it as the correlator of 
twist operators in CFT„. 

On the bulk side, to calculate the partition function on a higher genus Riemann surface 
holographically, one needs to find the gravity configuration whose asymptotic boundary is 
exactly the Riemann surface. To find the gravity configuration, one can use the Schottky 
uniformization to get the Riemann surface and then extend the uniformization to the bulk. 
However for one Riemann surface, there may be more than one Schottky uniformization and 
different uniformizations give different gravity configurations. Among all of the gravity con¬ 
figurations the one with the least classical action dominates the partition function in the large 
c limit |20| . The contributions from other configurations are suppressed as 0{e~^). Further¬ 
more the 1-loop correction can be determined by the functional determinant of the fluctuations 


around the classical background 26 by using the heat kernel method developed in 27,28 


Recently, by using the operator product expansion of the twist operators, the holographic 
computation of the double-interval Renyi entanglement entropy for the CFT has been checked 


beyond the classical level [^,29-^. Furthermore, for the single interval on a circle at finite 
temperature, if the interval is not very large, the holographic computation has been confirmed 
to be in exact agreement with the field theory computation , in which the thermal density 
matrix is expanded level by level [33| . 

In this paper, we study the Renyi entropy of a large interval on a circle at high temperature 


in the context of AdS 3 /CFT 2 correspondence, extending our previous study in 14 . In the 
large interval limit, the computations in both the field theory and the bulk need to be developed 


^In the literature, this n-copied CFT is usually called orbifold CFT. As in our following discussion we often 
use the tensor product of n copies of CFT, and we would like to call it the n-copied CFT. 
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furthermore. On the field theory side, the study in [14] showed that the perturbative series in 
the partition function do not converge well. Actually, the classical part of the Renyi entropy 
is just 


Sn {classical = ^ j—^ log sinh(27rry) + const. - ^ -?^<{ sinh'^(27rry)e 

bn 9 


+4sinh^(27rTy) cosh^(27rTy)e 
-ll + 28n2 + 199n^ 


fivrTR I “11 “ + 1309n'^ , , ^ . 

+ <--c„.h(i6,r,) 


w ^ s 77 - 346n2 + 197n'^ , , ^ , 

cosh(127rTy)- - -coshlSyrTT/) 

1440n4 ^ 2880n4 ^ 

-77 + 436n2 + 433n4 ^ ^ “77 + 466^2 + 907n^\ \ 

-cosh(47rry) +--J e 


1440n4 

^Q(g-107rTi?), 


2304n4 


( 1 . 10 ) 


where y is the length of the interval. When the length of the interval is comparable with the 
size of the circle y ~ R, the expansion converges very slowly and is not good anymore. This 
asks us to find another perturbative way to compute the partition function more effectively 
and reliably. In [^, we proposed to insert a complete twist sector states through the branch 


cut and expand the Renyi entropy with respect to e . In 


B 


15 


, we tested this proposal 


and reconsidered the noncompact and compact free scalars and found good agreements with 
direct expansions of the partition functions. Now we are going to consider the CFT with a 
holographic gravity dual, in the large interval and high temperature limit. We only consider 
the vacuum module of CFTi and its correspondents in the twist sector of CFT„. After cutting 
through the branch cut, the Riemann surface still have a nonzero genus and the four-point 
functions in the twist sector cannot be calculated directly. Furthermore, we use the OPE 
of the two twist operators at the branch points and compute the correlation functions on the 
unfolded cylinder of length n/I with the fields in the OPE at the different positions. We manage 
to expand the result with respect to complementary part of the interval length 1. We calculate 
the Renyi entropy up to order and e“^, including the leading linear c, the c-independent, 
and 1/c parts. Moreover, we find that the leading contribution is actually captured by the 
twist vacuum module. We support this result not only by the argument using the large central 
charge limit of the conformal blocks but also by direct computation using the Ward identity. 
As a result, we obtain the exact formula for the entanglement entropy. 


Eor the holographic calculation, we follow the treatments in p6] 20 , but basing on a 


different monodromy condition. As shown in 34 , the holographic entanglement entropy for 


the large interval case is not read from the bulk geodesic ending on the interval. Instead, 
it is the sum of the horizon length and the geodesic of the complementary interval. This 
fact suggests that there is a phase transition when the interval becomes large, and the bulk 
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Figure 2: Canonical cycles for the large interval at high temperature: n = 3 case 


gravitational conhguration for the large interval must be different. Instead of the cycles used 
in , we choose another n cycles to be of trivial monodromy. Among them, there is one cycle 
crossing the branch cut n times, and the other n — 1 independent ones crossing the branch 
cut and enclosing the complementary part of the branch cut in different sheets. As shown in 
Fig. i for n = 3, we set Aj’s to be of trivial monodromy. As a warm-up, we compute the 
classical holographic entanglement entropy by using the new monodromy conditions, and we 


obtain the result suggested in 34 . Furthermore, we compute the holographic Renyi entropy 
up to and for classical contribution, and up to and e ^ for 1-loop quantum 

contribution. The results are in perfect match with CFT’s computation. 

The remaining parts of the paper are organized as follows. In Sec. 2, we present the field 
theory computation. After a brief review on the twist sector of the CFT^, we focus on the 
vacuum module and compute the Renyi entropy in the first few orders. In Sec. 3, we show 
how to do holographic computation with the new monodromy condition. We obtain both the 
classical and 1-loop quantum results perturbatively. Up to the orders we are interested in, we 
find good agreements with the field theory results. In Sec. 4, we end with conclusions and 
discussions. We collect some technical details in Appendices. 
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(b) translate the interval 



Figure 3: Riemann surface for finite temperature Renyi entropy. The horizontal line denotes 
the spatial direction of unit length, and the vertical line denotes the thermal direction of length 
/?. (a) The Riemann surface for n cylinders connected by a branch cut, which is denoted by 
a solid line ■ In each sheet, there is a cycle along the imaginary time direction. 

However there could also be a cycle marked by which crosses the branch cut and goes 
from one sheet to the next one n times until it goes back to the original sheet, (b) The same 
as (a), just by shifting the branch cut to the boundary. denotes the complement of the 

original interval. In (c) we unfold the twist. There are n copies of cuts in the unfolded 

picture. 


2 Field theory calculation 


In this section, we present our computation on the large interval Renyi entropy on a circle at 
high temperature in the CFT which corresponds to pure AdS^ gravity. By the replica trick 
the Renyi entropy can be transformed into calculating the partition function on a higher genus 
Riemann surface, which is obtained by pasting n tori along a large branch cut. In Fig. 2, 
we cut open the torus and show the branch cut (interval). The horizontal line is the spatial 
direction of length R, and the vertical line is the thermal direction of a length (3. The large 


interval is presented as double solid lines between and ^ In Fig. 3b, we translate the 


interval, and in Fig. 3c, we unfold the branch cut and get a cylinder of thermal length n/3, 


with n cuts. We denote the coordinate in Fig. 3b as u and the one in Fig. 3c as y, and set 
the two branch points to be at in Fig. 


3b 


As in our previous paper [^, we would cut the Riemann surface through the branch cut, 


A^^'> cycle in Fig. 3b and insert complete bases in the twist sector of CFT„ at the cycle. The 
calculation transforms into a series of four-point functions with two twist operators and two 


9 











































operators in the twist sector. As the two branch points are actually very close to each other, 
we may take the OPE for the two twist operators, which amounts to an infinite summation of 
local operators in the n-copied GET. 

Alternatively, we may unfold the Riemann surface as in Fig. and insert complete state 
bases of the normal sector in single sheet GET to do the computation. Ghanging into the 


coordinate in Fig. 3c, we find that the localized operators in the jth copy in the OPE sit in 
the imaginary axis 

yO) = i(j/3). (2.1) 


In this way, each term in the expansion is a multipoint correlation function on an infinite 
cylinder, with two normal sector operators at the left and the right infinities of the cylindei|^ 
and n localized operators at . 

2nR 

We calculate the leading terms for the Renyi entropy expanded with respect to both e 
and 1. The results in the large c limit actually include the leading contribution, which is linear 
c, and the next-to-leading contribution, which is of order one, and even the next-to-next-to- 
leading contribution. We argue that only the descendants of the twist vacuum, the states 
generated by acting the Virasoro algebra Vir^*) of GFT„ on the twist vacuum, contribute to 
the leading c result. We confirm this fact by the analysis of the classical conformal block 
expansion in the large c limit. We furthermore derive the first few leading linear c terms 
from the descendants by using the Ward identity on the correlation functions of four twist 
operators, which could be related to the one in the double-interval case. In the re —)• 1 limit, 
we find that only the twist vacuum module has a nonvanishing contribution. This leads to 
the entanglement entropy of the large interval, which is exactly the same as the one from the 
holographic computation. However, when we consider the next-to-leading contribution to the 
Renyi entropy, we have to take the contributions from other states in the twist sector into 
account. 


2.1 Twist sector 

In this subsection, let us give a brief review on the twist sector in the GFT„. By the replica 
trick the Renyi entropy can be transformed into the partition function of a single copy GET 
on an re-sheeted surface connected at the branch cut. From the path integral, it is easy to see 
that the partition function could be taken as the re copies of field theory, one on each sheet, 
with fields on different sheets being related at the branch cut. In this re-copied theory, the 
locality requires us to introduce the twist field or antitwist field at the branch points [4 

^By the state-operator correspondence, the inserted normal sector states can be transformed into two vertex 

2-kR ^ 

operators at the left and the right infinities with a factor e "3 


35 
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( 1 ) 



Figure 4; Branch cut 


Let us show how the twist sector arises in a CFT„ with a branch cut, following the dis¬ 
cussion in [^. As the twist field is a local field, we consider simply the n-sheeted surface 
connected by a single branch cut. As in Figj^ the denote the branch point located at 
zi,Z 2 and the double line denotes the branch cut. Now we study it as a CFT,i. Considering 
an operator , it will change to (z) when it moves once around the branch point 

T'*' along the circle (2). The point z = zi is a branch singularity in the n-copied theory, which 
is a source of stress tensor 

<r(.)>=„<Ta)(.)>~A(„_i)^, p,2) 

when z is close to zi. Independence on the branch cut implies a local operator at the branch 
point. This local operator is known as the twist operator, denoted as 'T'^{zi). It is a primary 
field, with conformal dimension h = -^{n — In a similar way, we can get the antitwist 
operators at the other branch point. Much information of the CFT„ is encoded in the twist 
operators. For example, the partition function of an n-sheeted complex plane with N intervals 
is determined by the 2A^-point function of the twist and antitwist operators on a complex plane. 
Moreover, considering the operator-state correspondence, the twist operator corresponds to the 
ground state in the twisted sector of Hilbert space. Considering the OPE of the twist field 
with other basic fields in the theory, we will find other excited twist fields. Correspondingly, 
we find the excited states in the twist sector, as we will review soon. 

Before our discussion of the twist sector and antitwist sector states in the n-copied field 
theory, we show that the OPE of a twist sector operator and an antitwist sector operator give 
trivial sector operators in the n-copied field theory. As argued above, the excited states in 
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the twist sector could be obtained by considering the monodromy of the field moving around 
the branch point zi along the circle (2) in FigJ^ . However, when we consider the OPE of 
the operators in both the twist sector and antitwist sector, by the monodromy condition, the 
resulting states in the circle (1) in FigJ^ must be in the trivial sector. This fact has been 
applied in the discussion of the OPE of two twist operators in the short interval limit. In that 
case, the operators in the expansion are in the tensor product of the normal sector of the CFT 


in each sheet (which is just the trivial sector in the CFT„ picture), as shown in (12,25 . More 


generally, for the excited states in the twist sector, their OPE should consist of the trivial 
sector states. 


Now, let us give a review on the twist sector states 14 in the CFT„ arising from the 


replica trick in calculating the Renyi entropy. Let us work in the coordinate u in Fig. |3b[ 
To expand the partition function, we need to insert complete twist sector bases along the 
cycle We may temporarily forget about the geometric structure of the torus and only 

consider the geometry and the monodromy condition near the cycles . Moreover, as 

the vacuum module dominates in the large central charge limit, we focus on the twist sector 
from the vacuum module. For more complete discussion on other modules, please see for 
more details. In the vacuum module, the fields are constructed from the stress tensor. The 
monodromy condition on the stress tensor is 


(u + i/3) = (u), j = 0, • • • , n - 1 

in the u coordinate, with T^'^\u) = and 


(2.3) 


T{y + inl3) = T{y), 


(2.4) 


in the y coordinate. In the y coordinate, the n copies of fields are unfolded as 

r(^)(u)=T(y) (2.5) 


Taking the conformal transformations 


z = 

2ttu 

e P , 

(2.6) 


2-Ky 

(2.7) 

w = 

e "P , 

the monodromies in the new coordinates z, w 

are, respectively, 


r(i)(^e2^i) 

= r(^'+^)(z). 

(2.8) 

T{we^^^ 

) = nw). 

(2.9) 
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The states inserted at the cycle in Figs. 3b and Fig. 


3c 


can be described as the vertex 


operators being inserted at the origin of z and w. As in [^, we can redefine the operators in 
the z coordinate, 

n 

27vi I 


= fc = - ,n - 1 

i=i 


and expand it as 




r (^) 


m£Z ^ 


m+2+i 


( 2 . 10 ) 


( 2 . 11 ) 


(k) 

The operators Lm satisfy a commutation relation similar to the Virasoro algebra. Among the 
operators is of special importance. It is the total stress tensor for the whole 

n-copied theory, and {L^'\ are the generators of the corresponding Virasoro algebra Vir^*^ 
We may study the spectrum of the theory with respect to Vir^*). From the commutators 


between and 


[L®,LW] = -(m + ^)L^. 




n 


( 2 . 12 ) 


ik') 

we know that when the operators Lm act on a state, those with 

m > 0, 0 < k < n, or m = 0, 0 < k < n, (2.13) 

decrease conformal dimension, so they are annihilation operators; while those with 

m < 0, 0 < k < n, (2-14) 

increase the conformal dimension, so they are creation operators. Therefore we can define the 
vacuum for the twist sector to be 


^in I ^) = 0 for m > 0, 0 < k < n, or m = 0, 0 < A: < n. 


(2.15) 


The twist vacuum has the lowest conformal dimension 

f^v = ^n{l-^). (2.16) 

24 

Acting with the creation operators on the twist vacuum we can get all of the excited states in 
the twist sector. 

There is a one-to-one correspondence between the twist sector states in the CFT„ and the 
normal sector states in the original one-sheet CFT. Actually the trivial monodromy condition 
in the w coordinate suggests that the mode expansion in the w coordinate for the field gives 
the normal sector of the CFT. The conformal dimensions between the twist sector states and 
the normal sector states are related by 

1 cl 

htwist — -\- ^—71,(1 ^). (2.17) 

n 24 
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On the cylinder, the energy of the state could be written as 


H 


27r - (T f ^ 

+ ■^twist -I o ) — nK^'^ormal ^normal 

p 12 np 12 


(2.18) 


The ^ factor in the last equation is due to the fact that in the y coordinate the length of the 
thermal cycle is n/3. For convenience, we will denote the states in the twist sector as | t,i), 
corresponding to the state | i) in the original theory. 

It turns out to be more useful to classify the states in the twist sector by using the conformal 
symmetry in the re-copied theory. The states should be decomposed into different irreducible 
modules of the Virasoro algebra Vir^*^ To show this decomposition we calculate the chiral 
partition function for the twist sector 


Tvq^° = 

s=2 



• n 



) n 



(2.19) 


In the second equation, we have used the fact that there is a one-to-one correspondence between 
the twist sector states in the CFT„ and the normal sector states in the one-sheet CFT, with 


their the conformal dimensions being related by (2.17). For the vacuum module in the normal 
sector, the descendants are generated by the Virasoro algebra with s > 2, so the product 
begins from s = 2. For a primary operator with the conformal dimension h, if there is no null 
state in its descendants, its contribution to the chiral partition function is 


Trg^o 


1=9'' n 


1 

l-q^' 


( 2 . 20 ) 


Considering this fact, the quantity in the parentheses of (2.19) can be taken as a generating 
function for the primary operators with respect to the Virasoro algebra. As explained in [^, 
the operators in the z coordinate correspond to the generators T_mn in the unfolded w 
coordinate with mn > 1. As L-rnn do not generate null states in the normal sector vacuum 
module, the operators do not generate null states in the twist sector either. Therefore 
there is no null state in the descendants of each primary states in the twist sector. Expanding 
the function 


[q^^ • n 


1 




1 — 


-)i 


( 2 . 21 ) 


with respect to q, the coefficient before q^ is the number of the primary operator with conformal 
dimension h. It is clear that in the twist sector there are many new primary states and the 
number of the primary states increases exponentially with the conformal dimension. For 
example, acting with the operators on the twist vacuum, the resulting states that have 

conformal dimensions h = hy + ^ are the primary states, since they can be annihilated by the 
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operators , m > 0 in Vir^*). Among the modules in the twist sector, the vacuum module 
generated by , m < 0 on the twist vacuum is the most important one in our following 
discussion. We call this module the twist vacuum module. 

In the following discussion, we will meet another notion, the normal sector, in the n-copied 
field theory. It is dehned with respect to the n-copied field theory without a branch cut, or 
the tensor product of the n copies of Hilbert space of the normal sector of a single CFT. 
Moreover, as we only focus on the vacuum module of the CFT, we call the tensor product of 
n-copied vacuum module as the trivial sector in the n-copied field theory. Note that the stress 
tensor T^^o) jg g^jj^ well-defined quantity in the n-copied field theory without the branch cut. 
Therefore we may classify the states in the trivial sector by the Virasoro algebra Vir^*). In this 
case, we find that there are exponentially increasing primary operators with respect to Vir^*) 
in the trivial sector of the CFT„ as well. Considering the chiral partition function, we have 


OO ^ 

= <11 rr^)’ 


S = 2 
OO 


n 

s=2 


1-g' 


n 




- 1 


OO - 

S = 1 ^ 


( 2 . 22 ) 


In the last equation, we decompose the whole partition function into the contribution from 
different modules with respect to Vir^*^ Each module is generated by > 0 acting 

on the highest weight state. The hrst term denotes the module generated from the vacuum 
state, with zero conformal dimension, so the product starts from s = 2. For the other primary 
operators, there are no null states in their descendants. This is because the primary state 
has nonzero conformal dimension, and the states \ h) have nonzero norm. Considering 
[Lr,L-r\ = ~ 1) + 2rLo, for r > 1, the commutator has a linear c term. In the large 

c limit, all of the states have nonzero norm. To read the number of other primary states, we 
just need to expand the quantity 


n 


1 

1 


-1 


(2.23) 


vii (1 - / 

with respect to g, such that the coefficient before g^ is just the number of the primary states 
with a conformal dimension h. We can easily see that the number of the primary states 
increases exponentially with their conformal weights. In the following, we call the module 
generated from the vacuum a trivial vacuum module, with respect to Vir^*). 


2.2 Renyi entropy 


As discussed previously, we can expand the partition function by inserting complete twist 


sector bases in Fig. 3b or complete normal sector bases in Fig. 3c The partition function can 
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be expanded as 


I r-(ni)r+(n2) I 

i 

= I f~iui)f^{u2) I t,i) |„, 


(2.24) 


in the u coordinate, and in the second line we use (2.17). And the thermal partition function 
reads 


= 


E' 


'n = / e 

i 




(2.25) 


In the large c limit, we only need to consider the vacuum module, which captures the 
perturbative effect. We list the first few states in the vacuum module and their vertex operators 
at the origin and the infinity in Appendix A. Such states give the first few leading order 
contributions to the Renyi entropy. 


Each term in (2.24) is a four-point function with two twist operators at the branch points 
and two operators in the twist sector at the left and the right infinities. For the CFT 
dual to pure gravity, we do not know exactly the analytic form of this correlation function. 
Nevertheless, when the twist operators are very close, we may take the OPE of two twist 
operators 




(2.26) 


where 


viu) = i + ^l(i-l)^2rW(u) +j;^(i-l)zW*)(n) 

12 488 

i i 

+ E ■■ -^aV‘\u)) 


+ Ei 

^ l<^2 


1 1 


1 


8cn4siiP^(n-ia) 144 


p- + J^(l - + O{l0), (2.27) 


and the similar form for the antiholomorphic part V{u). In the operator product expansion, 
we only consider the vacuum module in each sheet and ignore the other modules. 

In practice, it is more convenient to unfold the twist and consider the correlation function 


in Fig. 3c After transforming into the coordinate y, we find that the operators in R(0) are 
localized at 

rO)(0)^r(y)|,=,,-^, (2.28) 


■^The bases we insert when cutting the Riemann surface are in the Schrodinger picture. In (2.241, we change 

27rR jj - 

them into the Heisenberg picture, with a factor e a , and the states now correspond to the vertex operators 
at the infinities. 
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namely, the operators in different sheets are unfolded and located at different positions in the 
cylinder. We may use the OPE of the twist operators to compute the partition function per- 
turbatively. Now the partition function can be expanded as a multipoint correlation function 
on the cylinder involving n operators located aX y = ijj3, and two vertex operators at left and 
right infinities. We can furthermore take a conformal transformation into the w coordinate, 
and do the calculation in a full complex plane. Formally, we still write 


77 Rc A. ^ ^ 

Zn = {t I T~{ui)T^{u2) I t)(y^ 


I f (ui)f'+(u2) I 

{t\f-iui)f+{u2)\t) 


C /-l _ 1 \ 

= Cnl ^ ^ ^ 


{t\v(o) \t){t\v(o) \t)(y ^ \ \ : \ I e 0 (> ' ) : \ I e 

' ' ^ ^ (t I y(o) I t) {t I ^(0) I t) ’ 


= e^^Pcnl E(0) 


ttRc c^/i 1 \ 


(2.29) 


In the second line we take into (2.26), and the fact that the holomorphic and antiholomorphic 
parts decompose. In this case, the holomorphic and antiholomorphic parts are equal to each 
other. The Renyi entropy is 


Sr,. = - 


1 , Zn 

log 


n — 1 ° Z” 

1 c _ In,, ttRc 

--((logcn - -n 1- 2 log^ + 

n — 1 6 bnp 


+ 21og(t I R(0) I t) 


, (^5 * I 1^(0) I T *) -^/i NN .TtRc , ,v-^ _22 lR;, 

^ '))-«(-^ + 21og(^e /3 '))). (2.30) 


Denote 


fi = 


{t,i I P(0) I t,i) 

(t,0| R(0) |t,0) 


{i I P(0) I i) 


(0 I P(0) I 0) 


(2.31) 


where in the last equation, we change into the y coordinate. Taking /j into (2.30), we get the 
first few terms of the nth Renyi entropy 

2 


Sn = - 


1 , c n + 1 , , ttLc n + 1 

-logc„ + --log! + —- 

“ 6p n 


n — 1 


log(t I R(0) I t) 


n — 1 6 n 

2 4-7, fl e-n-R ^ _ __ __ 

--(/ 2 e + fse + (/4,1 +/4,2 --/ 2 )e +0{e )) 

n — 1 1 

2n 


X 2 _ StvR 


_ 477 J? _ 677fl 3 877fi _ 1077 J? 

+- -{e +e P +-e p + 0{e p )). 

n — 1 1 


(2.32) 


_ 277fi _ 277R 

The Renyi entropy is expanded with respect to e P and e . The two expansion argu¬ 
ments are not independent. In the above formula, we have actually done the computation in 


Fig. 3c In other words, we have unfolded the twist and consider the insertion of the normal 
sector states at the left and the right infinities of the cylinder. Meanwhile to calculate the 
analytic form of /j, we also take the OPE of the twist operators which is an expansion with 
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respect to the relative length of the two twist operators Therefore, we have two kinds of ex¬ 
pansion, one from the normal sector states and the other from the OPE of the twist operators. 
This results in two different expansions in Sn- 

The explicit expressions of /j’s can be found in Appendix C. There are a few remarkable 
properties on /j’s. 

1. First of all, the P terms in /j’s always take the form i ■ This is because 

the operator at P order in the OPE is the stress tensor T{z) = whose 

correlation function is fixed by the Ward identity. 


2. Second, for each fj there exist some exceptional integers re at which fj does not share 
a general formula and take specihc form. This fact forbids an analytic continuation of 
the Renyi entropy to noninteger re in order to calculate the entanglement entropy. It is 


because in the OPE (2.26) there is always a term like 

jl <32 

The correlation function involving such a term includes the summation 

n—1 


(2.33) 


E 

ii-j2=i 


7r(/i -/2). 


re 


(2.34) 


When m < 0, the summation has a universal formula for any integer re, while for m > 0, 
it is more complicated. We can rewrite the sine function in the summations as 


. h)-. sr^ 2-^ir(n-n) 

sm (-) = > are 

n 


If we take a summation for each term we find 

n—1 


E » 

ii-i2=i 


2-!rir(ji-j2) 


-1, for ^ ^ N 

re — 1, for ^ £ N ’ 


(2.35) 


(2.36) 


which shows the nonanalytic origin. 


3. We notice that fj has no linear c contribution for nearly all of re but finite exceptions. 
Actually, when | t,j) is primary, there is no linear c term in fj. To understand this 
effect, we can first transform into a z coordinate and expand the four point function by 
conformal blocks. Using the large c conformal block [^, there should be no linear c 
contribution. For nearly all of re but finite exceptions, | t,j) is primary. For example, 
when re is bigger than the conformal dimension of | j), which is the corresponding state 
in the normal sector, the state | t,j) is primary because it can be annihilated by all of 
the Virasoro algebra generators with m > 0. 
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4. Furthermore, we also notice that for the n in which fj has no linear c terms, fj always 
share a general formula up to As we discussed before, the non-analytic property 
comes from correlation functions 

(2.37) 


Because of the symmetry, the correlation function is a sum of 

sin--(ji-j2). (2.38) 

n 

If there are some terms with m > 0 the final result is non-analytic for n. However, if 
I t, j) is primary, this cannot happen. Consider 


^ {t,j I I t,j) 

3i<32 

n—1 

= ^ X] I I t, j), 

ii-j2=i 

n—1 

= \ Y. J^femSin”*^(jl -J2). 

ii-i2=l rn 

On the other hand, we have 


h<32 


- n—1 - n—1 

2(j;rO)(u))2-i^(rW(u))2 

j=0 j=0 





(2.39) 


(2.40) 


For the primary state, the correlation function {t,j \ T{u)‘^ \ t,j) can be fixed by the 
Ward identity, and it is analytic for all n. The term {t,j \ {T^^\u))‘^ \ t,j) does not 


depend on j, and it is also analytic, which means there is no m > 0 terms in (2.39). It 
is not clear whether this property can be extended to a higher order of OPE expansion 
with respect to 1. 


Here we just list the first few leading order results of the Renyi entropy: 


( li , .TT^P 




■)- 




64/3^ 


+ - 


+( 1 - 


2048/34 

37r2/2 


“ (2 + wtitv + 


2/32 ' 384/3' 

2757r4/4^ 


184320/34 

257r4/4 _ 2 ^ 37r2/2 

-^)e -(2 + ^7^ + 


ll7r4/4 

e 


4/32 256/3' 


157r4/4 3nR 

- ^ 


2/32 


384/34 


-)e~ 


1 


+ - - 


_2zfl 7r4/4 _3^ 2l7r4/4 

^3 + 


64/34 


128/34 


64/34 


+0(/^) + 0(e"~) 


(2.41) 
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^3 


1 , 27rR ^ 3A7 ttH^ , 128tt^I^_mr\ 

“ 2 ^^ ^ ^ Vo ^ ~W ^ 243/32 “ 1771470^4^ “ 6561/34 ^ 

f levr^Z^ . IGtt^I^ SOOtt^/^, .UlR , Stt^I^ 2807r^^^, 

1 177147/34 ^ ^ 81^2 ^ 59049^4^ V ^ 27/32 6561/34 ^ 

,3 167r2/2 6087r4/4 _ 82 rfi'| 1 / 2567r4/4 327r4/4 

^2 27/32 6561/34 ^ i cV 59049^4 6561/34 

- IOttR 

+0(l^) + 0{e-^) 


_27rR 

e ^ 

647r4/4 stvr \ 

-e 3/3 

6561/34 y 

(2.42) 


54 


1 ^ 1 / 20897r4/4 1267r4/4 _2^\ 

--C4 + C log + + ( 2304^2 “ 14155776/34^ “ 6144/34 ^ 

f ll7r4/4 2 57r2/2 7497r4/4 2 bir'^P 13757r4/4 _3^ 

\ 262144/34 ^ 3 487r2 147456/34^ ^3 32/32 98304/34 ^ 

, 57r2/2 7197r4/4 _22rfi'| l / 197r4/4 _2rH 957r4/4 237r4/4 _22rfl\ 

^ I 6/32 16384/34'^ / c V 24576/34 49152/34 8192^4 y 

, b'KR 

+0{f) + 0{e~^) (2.43) 


5. 


n — 1 


log Cn + c - 


i n -I- i 


6 n 


log / + 


/I J t /t- T 


6/3 n 


+ 


36/32n3 


/47r4(n + l)(n2 + ll)(3n4 + lOn^ + 227) 
226800/34n7 


n — 1 


12960/34n7 


1 + 


27r^/^(n^ — 1) 


3/32n3 


TT 


+ 1 + 


/4 (n® - 12n3 - 3n4 + 24n3 + 135n2 - 444n + 299) \ _ izLfl 

---- P nP 

54/34n'^ ) 

72 (n2 - 1) 7 r 4 / 4 (n 6 - 18n3 - 7>v* + 36n^ + 835n2 - 4626n + 6081) \ 

B'^n^ 36/34n’^ j ^ 


+ - 

for n > 4. 


/32n3 36/34n’^ y 

r2/2(n2-l) /47r4(5n‘^-120n3-31n4 + 240n3 + 9635n2- 72840n + 135111)\ _ 82 l «1 

---- — ---- \ P "/5 ; 

/32n3 90/34n'^ / j 

/^7r4/4(2n6-63n4 - 672n2+ 733) 7r4(2n6 - 217n4 + 1148n2 - 933) 

I -e "/^ H-e "/3 

o \ 945/34cn’^ 630/34cn'^ 

/47r4(2n®-427n4 + 13328n2-40320n +27417) _io 2 l«, 

^-'e «)+0(/‘) + 0{e 


+ (2 + 

1 2 
+-; 


(2.44 


2.3 Classical limit of the conformal blocks 

In the previous subsection, we claimed that for each fj, when | t,j) is a primary state in 
n-copied theory, it has no linear c contribution. In this subsection, we clarify this fact from 
the point of view of a large c conformal block, and furthermore we show that we only need to 
consider the twist vacuum module to find the linear c order entanglement and Renyi entropies. 
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Let us study the four point function between two twist operators and two vertex operators 
corresponding to | t,j) in the twist sector, which is primary under the conformal symmetry of 
the n-copied theory. By the conformal transformation 

v = eP^ ( 2 . 45 ) 


the fonr-point correlation function can be transformed into 

{t,j 17'“(-^)r+(^) I t,j) \u 

27r 27ri c / 1 \ 27ri , 

= (ye /3 (e / 3 )r+(l)|Lj) 

which can be expanded by the conformal blocks as 


E 


(jhp 


2tiI 

F{nc, hr, hr, htj, hj] hp-,1- e p ), 


( 2 . 46 ) 


where nc is the central charge of the CFT„, and CjJ'® is the OPE coefficient from two primary 
operators with hi /i 2 conformal dimension to a primary operator with /is conformal dimension. 
The first four conformal dimensions in T are for the four external operators, two twist operators 
and two operators in the twisted vacunm, and the last one hp is the conformal dimension of 
the primary field in the propagator. In each replica, we consider only the vacunm module in 
CFTi, so the states in the propagator are in the tensor product of n vacuum modules, which 
is the trivial sector in the CFT„. 

One essential point is that the OPE coefficient is of order . The primary operators can 
be normalized as 

(0(zi)0(z2)) = (2-47) 

In our case, each operator in the propagator is a combination of the stress tensors and their 
partial derivatives. If the largest number of the stress tensors in the combination is r, snch 

r 

an operator should be normalized by a factor of order c “2 in the large c limit. The OPE 
coefficient equals the expectation valne in the n-sheeted surface. To compute it we need 
to transform into a full complex plane w. We can decompose the transformation into two 
steps: the first one transforms z coordinate into an n-sheeted fan with boundary condition 




( 2 . 48 ) 


and the second step nnfolds the n-sheeted fan into the full complex plane. In the transforma¬ 
tions, the number of the stress tensors in the operators does not change, so the expectation 
valne is at most order in the large c limit, which means the OPE coefficient is at most 
order 


25 


. Eor Ojj /j, we jnst need to insert two extra operators, and it is still of order (^. 
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Furthermore, the leading contribution in the conformal block T is the same in the large c 


limit. As suggested in 21,37 


T{c,hi,hp,x) =exp[-^/(/if\/i(,°),x) + 0(c°)], 

where 

= hm ^ 

= lim —. 

^ C-^OQ C 

In the case at hand, we have the relation 

hp = 0(c°) ht,i -ht = 0(c°), 


(2.49) 


(2.50) 


(2.51) 


so that the classical conformal blocks are the same for all different terms in the expansion. 


Taking into (2.49), it is easy to prove that the four-point functions are independent of j in 
the leading order, if | t,j) is a primary operator. Moreover, even for the two twist operators’ 
correlation on torus, one can also only consider the twist module generated from | t,j) by 
Vir^*) for the linear c order. Other modules gives only ^ corrections. 


2.4 Leading contribution from the twist vacuum module 


As we showed above, it is only necessary to insert the twist vacuum module to compute the 
leading c order Renyi entropy. In this subsection, we use the Ward identity to calculate the 
contribution from these terms explicitly. By using the Ward identity, all of the multicorrelation 
functions for the descendant operators can be derived from that for the primary operators. 
From the recursion relation 


i 

m 

= E<oi( 


t=i 


2T{wk) dT{wk) 

~r / \ o 


k=2 


2{wi - Wk)^ (rci - Wk^ wi-Wk' 






t=l 


n m 


, f hk4>^’‘\zk) , d4>^^\zk)~^ y-I- (j)^ 

+E(o 111 


k=l s=l 


t=l'^n,t^k 


where O in the bra denotes a primary operator in the infinity. With proper contour integral 
and contraction, we can derive any correlation function of the descendants of the primary 
operators. 

What we need to compute in the partition function are the ratios 


{t,i' \T {ui)T~^{u 2) \t,i) _ {t,i \T {zi)T~^{z 2) \ t,i') 

{t\T-{ui)T+{u 2 )\t) ~ {t\T-{zi)T+{z 2 )\t) 
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Here | t, i') means the Virasoro descendants in the twist vacuum module generated by acting 
plane by 


> 0 on the twist vacuum. And we have changed the coordinate into a full complex 


21, 

z = ei^ 


(2.54) 


so that the two inserting operators are at the origin and the infinity respectively. By the 


conformal transformation and the Ward identity, all of the terms in (2.53) can be calculated by 
the four-point functions, which are related to the double interval mutual information. Actually, 
for the simplest case, the contribution of the twist vacuum is encoded in the correlation function 
of four twist operators with two of them being inserted at the origin and the infinity in the 
complex z plane. In general, the correlation functions of four twist operators read 

log < r+(zi)r-(z2)r+(z3)r-(z4) > 

= (n - + log < r+(zi)r-(z4) > +log < T-{z 2 )T+{z 2 ) > 


= (n - - :^{n - -)log((zi - Z4)(^l - Z4){z2 - Z3){Z2 - Zs)), 

12 n 


(2.55) 


where 


X = 


^^3 - 2^2 ^4 - Zi 


Z3 - Zl Z4 - Z2 

and is the mutual Renyi information. If we set one point to infinity, then we have 


(2.56) 


log<t|r+(zi)T-(z2)r+(z3)> 

= log lim 

24 —^CXD 


< r+(zi)r-(z2)r+(z3)r-(z4) > 


= {n-l)I^'^\x,x) - ^{n- ;|()log((z2 - Z3)(Z2 - Z3)), 


where 


x = 


n 


Z 3 - Z 2 


(2.57) 


23 - zi 

The perturbative computation of has been done in 


25 


(2.58) 

29|. We list them in Appendix D. 


With these results, we can derive any /■> for the descendants of the twist vacuum. 
The first few lowest descendants in the twist vacuum module are 


where 


a = — 


m = 


\t,l) 

r ^ 1 r(0) 1 

= (— (n -)) 2 L 

t) 



= 

1 i) 


t,2,2) 

= m(LL°jL® + oLL”^) 

1 t) 

(2.59) 

3(n^ 

-1) 




8n2-2 ’ 

c^(n^ — 1)^ 
72n2 


+ o 


[n? — 1)(5 -|- 7n^) 
24n(4n^ — 1) 


(2.60) 
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Note that these states consist of a special set of excited states in the twist sector, which are 
the decendants of twist vacuum generated by And their contributions to the leading 

linear c order are, respectively. 


/i' = 


f2,l' ~ 


( 2 _ ^^3 ( ,4 , (-11 - 2 n ^ + 19 n ^) 7 r ' ,6 

^ ^ Vl08/3^n5 ^ 9720/36n9 

(f,2,r |r-(-^)r+(^)|t,2,h) 

{t\r-{-!,)T+a)\t) 


4^_6 \ i^2 

l^ + 0{f)) +(1 + 


7\\ , n 1)’’'^,2 , 


2 2 + o{n) + oic-^) 




= c (—1 + n"^) 


2 n 2 A-1 + 4n^)^^,4 , (ll-9n^-llln^ + 289n6)7r® 6 7 , On 


54/3^n® 


r + 


4860/36n9 


-r + 0{V) +0(c9) 


/2,2' ~ 


= C 


(t,2,2Mr-(-|)r+(^) |t,2,20 
(t|r-(-^)r+(^)|t) 
,7(n2-l)V 


'23328/38nio 


2,^3 TT^ ^4 (—11 — 32n^ + 49n^)7r® 

^54/3^n^ 4860/3®n9 


+ 0(1^)) + c(-l + n^)3 ( ^ r: JJ ^ + Oif) ] + 0(C°). 

(2.61) 


The leading c contribution for the Renyi entropy reads 




(class) 


ttcR n + 1 n + 1 


+ C— — log l + c{n + 1) 


TT 


2 TT^in^ + 9n^ - 9n^ + 11) ,4 


rr - 


6/3 n ' “ 6n y36n^f3^ 12960n7/34 

7r®(19nio + 85n3 - 125n6 + 251n^ - 274n2 + 188) 

2449440nii/36 

/ .t^ 4 /'ll 0.^2 

+c(-l + n)^(l + n)3 ( - 


+ 


3/ ^4 (-ll- 2 n^ + 19 n^) 7 r 6 ^g\ _2^ 


54/3^n® 


4860/36n9 


+c(-l + n)(l + n)^ ( - 

67rR 

+0{f)+0{e- — ), 


{3v? - l) 7 r 4 4 (11 + - lOln^ + 209n®)7r® . 


18/3^n^ 


-r- 


1620/36n9 


^ \ 4l:7r/l 


(2.62) 


and 


^2;2 


„(class) _ , r-class _ ^ 

n^i ^ 3/3 Vl8/3^ 540/34 ' 8505/34 


7r4/4 71-4^6 

+ ) +0(/3), 


(2.63) 


for leading order of the expansion. 

From the result, we find that in the entanglement entropy there is no finite size correction 

2ttR 

proportional to the powers of e ^ . Such a correction, if it existed, should come from the 
4-point functions of the descendants 


(t,/,i|T {ui)T"^{u 2) \u= {t,i,i\T {zi)T^{z2)\t,i,i) |^, 


(2.64) 
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where | t, i, i) is a state generated by a set of creation generators acting on the twist 

vacuum | t). Consider the Ward identity 


[Lr,T {ui,ui)] = {r + l)hru[T {ui,ui) + u'l'^^dT 


(2.65) 


where 


hT= -)■ 

24 n 


( 2 . 66 ) 


In the right side of (2.65), because h-j- = -^{n — the first term is of order n — 1; because 


the four point function is constant when n = 1, the second term should also be of order n — 1. 
Similarly when moving L_r to the left side, the commutation term will contribute an n — 1. To 
calculate the correlation function, we move all of the annihilation operators to the right side 
and then move the reduced creation operators to the left side. It turns out that the leading 
contribution terms are at least of order (n — 1)^. Therefore 

c.(class) ^_ 

n - 1 1 6n/3 


+ log(t I T iui)T~^{u 2 ) I t) U ^ + 0{n - 1) 
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6j3 

17i^(^ - n) + log ^(^)f(^-^)(t I T~{zi)T^{z 2 ) I t) 1^^ I + 0{n - 1) 


n — 1 ( 6/3 n 
(n + l)7ri?c 1 


6n/3 

(n + l)7ri?c (n + l)c 


n-l 




|i(n—-) 

— h ' n ' 


+ 0{n — 1)^ > + 0(n — 1) 


Zl - Z2 


6n/3 


+ 


111 


6n 


logsinh — + 0{n — 1). 


(2.67) 


In the third equation we use the classical conformal block 21 , and the entanglement entropy 
is 

'TrnT-^ n tt] 

( 2 . 68 ) 


^(class) _ ^ c , ^ 

i=EE - 3 ^ +3 logsinh^. 


which matches with previous result (2.63) up to order This is the high temperature entan¬ 


glement entropy for a large interval and it satisfies the relation 


lim See{R - e) - See{() = Sth- (2.69) 

3 Holographic Renyi entropy 

In this section, let us calculate the entanglement entropy and the Renyi entropy holograph¬ 
ically up to 1-loop order. In the field theory, by the replica trick the Renyi entropy can be 
transformed into the partition function on a higher genus Riemann surface. Holographically, 
this partition function can be computed in the semi-classical AdS gravity in the large central 
charge limit. Based on the AdS/CFT correspondence the gravity configurations must be the 
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classical solutions with the asymptotically boundary being the Riemann surface [38] . More¬ 
over, for the same Riemann surface, there may be more than one gravitational solution. The 
partition function is the summation of the classical contributions and the quantum corrections 
at different saddle points. Among different saddle points, the one with the smallest action 
dominates the contribution, and other saddle points give non-perturbative corrections of order 
e~'^. Therefore, in the large c limit we only need to consider the saddle point with the small¬ 
est action. The regulated on-shell action of this saddle point gives the classical contribution, 
corresponding to the leading linear c result in the field theory, while the 1-loop determinant 
of the fluctuations around the saddle point gives the quantum correction, which corresponds 
to the order results in the field theory. 


As in 20 , we assume that only the handle-body solutions contribute to the partition func¬ 


tion. The handle-body solutions could be obtained by extending the Schottky uniformization 
of the Riemann surface to the bulk. In this section we first give a brief review on the Schottky 
uniformization and the on-shell action. Then we discuss the monodromy condition for the 
n-sheeted torus pasted along a single large interval to find the uniformization. We compute 
the classical part of the holographic Renyi entropy perturbatively. Furthermore, after carefully 
studying the primitive class of the Schottky group, we calculate the 1-loop corrections to the 
entropies, following the treatment in [28] . 

3.1 Schottky uniformization and the partition function 

In three dimensional AdS pure gravity, all solutions with constant negative cnrvature are 
quotients of the AdS space. In terms of the Poincare coordinates, the AdS space could be 
described as an upper-half space with the metric, 

dudu + d^^ 


ds^ = 


C >0, 


(3.1) 


where u = x + iy is the coordinate of a complex plane. The isometry group of AdSa is 
PSL{2,C) [^. The coordinates u,^ can be combined into a quaternion = u + j^, on which 
the isometry group acting as 


c = 


aC + b 


(3.2) 


cC + d 

with ad — be = l,a,b,c,d £ C. At the asymptotic boundary ^ 0, the transformation is just 

a linear Mobius transformation on a complex plane, 

au + b 


u 


cu + b 


(3.3) 


Generally, the gravity solution can be written as AdS’g/T, where AdS'^ is the global AdS^ with 
some fixed points being removed and T is the discrete subgroup of PSL{2, C). The asymptotic 
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boundary is r^/r, where 0 is a full complex plane with some fixed points being removed. If 
we focus on the handle-body solutions, the subgroup is just the Schottky group. 

Every compact Riemann surface can be realized by the Schottky uniformization. For a 
genus-sf Riemann surface X, its fundamental group is generated by 2g generators, 

ai, 02 , ... Og; 6 i, 62 , bg (3.4) 

with constraints 

9 

JJ[ai,6i] = l, {a,h\ = aha~^h~^. (3-5) 

One can always choose g loxodromic generators Li, L 2 , ■■■ Lg and a fundamental regain D 
bounded by 2g circles Ci, C 2 , ... Cg and Ci, C 2 , ... Cg, such that Li{Ci) = C^. Identifying g 
pairs of circles by the generators, we obtain a quotient space, which is just a genus-g Riemann 
surface. Here Oj is just the image of Ci{C[) under the quotient map in the homology group; 
the group of covering 0 —)• fl/E is the smallest normal subgroup J\f containing the elements 
Oj’s, and the Schottky group is isomorphic to 'Ki{X)/N. The Schottky uniformization can be 
extended to the bulk, which is an automorphism of the AdS^ space, with the Oj cycles in the 
bulk being contractable. 

For one Riemann surface, there are more than one way to choose the a* and bi cycles. 
Different choices of the generators of the fundamental group correspond to different realizations 
of the Schottky uniformization. Even though different Schottky uniformizations describe the 
same Riemann surface, their extensions to the bulk give different gravity solutions. 

The Schottky uniformization problem for a general Riemann surface could be solved by 
considering the differential equation 

= 0 (3.6) 

where R^^\u) is the Schottky projective connection on a marked Riemann surface. R^^\u) 
is uniquely determined by the normal subgroup M. Namely it depends on the choice of the 
generators. A ratio of the linearly independent solutions of the above equation determines 
the quotient map in the covering space D. More importantly, it turns out that up to a 
normalization R^^'^ is just the holomorphic stress tensor T{u) of the Liouville theory 
which is the regulated on-shell action of the bulk solution of the AdSs gravity. The explicit 
forms of the stress tensor depend on {3g — 3) complex accessory parameters with respect to 
the holomorphic quadratic differentials on the Riemann surface such that the determination of 
the uniformization map is usually a very difficult problem. However, for the Riemann surface 
in computing the Renyi entropy, the uniformization problem could be solved perturbatively in 
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some cases due to the replica symmetry. For the double interval case 20 , the stress tensor 
takes the form 


r(S)(u) = Y, 


A 


+ 


7* 


{u — Ui^ U — Ui' 


where 


A = 7 ( 1 - 4 )’ 


(3.7) 


(3.8) 


2 ' 

and there is only one conformal invariant accessory parameter. For the single interval on a 


torus, the stress tensor takes the form 26 


R^^\u) = ^(Ap(u - Ui) + 'jiCiu - Ui)) + 6, 


(3.9) 


where p is the doubly periodic Weierstrass function 


1 


p{u;R,-) = E I 

(m.n)^( 0 , 0 ) 

7r2r2 


1 


1 


E 


(u + mi?+fl)2 {mR+^^y 


^ ^ CIT 


sinh^ [ttT {u + mR)] sinh^ (vrmT R) 


+ 


VT' 


'.'-p2 


and 


Q{u) = ttT coth[7rr(ti + mR)] + — 

• ^ ^ Cl 








sinh^ nmTR 


(3.10) 

(3.11) 

(3.12) 


To solve the problem, one has to impose the monodromy condition on some cycles to fix the 
accessory parameters. The different choices on the cycles with trivial monodromy give different 
Schottky uniformization. 

On the other hand, the regulated on-shell action of a AdSa gravity solution is the so-called 


Takhtajan-Zograf action 38 . Moreover, the dependence of the action on the moduli parameter 


has been studied in 39,401. For the gravitational conhguration dual to the n-sheeted Riemann 


surface, the action obeys the equation 20 

dSn 


cn 


dui 6(n —1)"*"*' 

This equation allows us to obtain the classical action of the gravity solution corresponding to 
a Schottky uniformization. Among different uniformizations for the same Riemann surface, 
the one leading to the least gravitational action dominates the partition function. 

Here let us focus on the case that there is a single interval on a torus. Because the cycles 
around two branch points are always of trivial monodromy, we have 


(3.13) 


71 = - 72 - 


(3.14) 
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For convenience, we redefine the functions and rewrite the stress tensor as 


R^^\u) = Ap(u - ui) + Ap{u - U 2 ) + 'yC{u - ui) - 'yC{u - U 2 ) + 6, (3.15) 


where 


Piu) = ^ 




^ sinh^ [ttT (u + mR)] ’ 

00 

c(-) = E vrT coth[7rT(M + mR)]. 


(3.16) 


For the classical partition function, we need to calculate the on-shell action of the gravity 


solution with proper boundary terms as regulators. It turns out Eq. (3.13) is not enough 


to determine the action completely. Besides the dependence of the action on the accessory 
parameter, we have to take into account its dependence on the size of the torus. In [l4|, we 


proposed another differential relation on the partition function, in addition to (3.13), in order 
to determine the size dependence of the partition function completely 

dSn 


OR 


c n 
127r n — 1 


PiS-6n=l). 


(3.17) 


With Eqs. (3.13) and (3.17), we can determine the partition function completely. 


From the holographic entanglement entropy of one single interval in the black hole back¬ 


ground 34 , there should be a phase transition when the interval becomes large enough. This 


means that for a very large interval one should impose a different set of monodromy conditions, 
which leads to different Schottky uniformization. To support our choice on the monodromy 
conditions for the large interval case, we will compute in the following the holographic entan¬ 
glement entropy and compare it with the result in |34| . 

Let us first review the holographic computation in the short interval case. We set the branch 
cut at [ui,M 2 ]- At a high temperature, the thermal cycle should be of trivial monodromy, and 
the wave function transforms as 

ip{u+^) =--i/jiu). (3.18) 

If we transform into the z coordinate, there is no minus sign in the monodromy condition. As 


discussed in 26 , to compute the holographic entanglement entropy we only need to study the 


solution near n = 1 and expand the wave function and the parameters with respect to n — 1 
as 

00 

ip{u) = V'*'°^(u) -|- ^(n — 1)V*'*^('*^)) (3.19) 


i=l 


7 = 

i=l 




(3.20) 
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with 


5 = -27r2r2 + ^(n- lyP, 

i=l 


(3.21) 


Expanding the trivial monodromy condition with respect to (n — 1), we have 


at each order. With proper redefinition of A and B in Eq. (3.22), we can also set 
,/; 0 )( 0 ) = = 0 , for j = 1 , 2 ,... 


(3.22) 


(3.23) 


(3.24) 


Taking the expansions of the wave function and the parameters into the equation, we hnd the 
following equation at the leading order 


(u) — 'k‘^T‘^'4^^^\u) + -m(u)V'^°^(ri) = 0 , 


where 


m 


l{u) = p{u — Ml) + p{u — M 2 ) + 7^^^C(m — Ml) — — U 2 ) + 5^^^. 

We get the solution 

_ 'UjTT'J. ' /*Uj 

= -—— / e^^'^m{r)'ip^^\r)dr — -—— / e~'^'^'^m{r)Tp^^\r)dr. 


2ttT 


/o 


2ttT 


'0 


Furthermore, considering the last two equations in (3.24), we get 

rl 


/o 


IS 

m{—)ds = 0 , 




/o 


e^^^^mC-^)ds = 0 . 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


Solving these equations, we find 


= 0 , 

7(^1 = 27rTcoth7rr(M2 — Ml). 


(3.29) 


Taking 5 and 7 into Eqs. ( 3.13|3. W ), we obtain the classical entanglement entropy 

= g logsinh^ ttT{u 2 — ui) + const. 


(3.30) 


which is the geodesic length in the bulk connecting two branch points. 

On the other hand, for the large interval case, we choose another n cycles of trivial mon¬ 
odromy. To compare with the small interval case, we set the branch cut at [0, mi] |J[m 2 , 7?]. 
Among n trivial cycles, there is one cycle that goes across the branch cut for n times. This 
cycle is denoted as A^^^ in Fig. 2 and Figs. There are other n — 1 independent cycles 
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enclosing [ui,U 2 ], the complementary part of the large interval. These cycles are denoted as 
i > 1 in Fig. As in the small interval case, we expand the wave function and parameter 
with respect to (n — 1 ), 


■(/; = 'i/)(°)(tt) + ^(n — 

(3.31) 

2=1 

(3.32) 

00 

7 = 

(3.33) 


i=l 

and the zeroth order wave function is 

= A + B e~~^ 

The monodromy condition for the cycle is 

777 

'tp{u + —) = -tpiu). 


(3.34) 


(3.35) 


where u + y means that the argument goes across the branch cut for n times. Expanding 
the wave function with respect to n — 1, we find the same differential equation (3.25) at the 
leading order, and the same and 7 ^^^ as in the short interval case. The only difference 
comes from ■ Taking them into Eqs. (3.13) and (3.17), we obtain the classical 

entanglement entropy of the large interval 


^^chissical) ^ c 


(3.36) 


Namely, the holographic entanglement entropy(HEE) of a very large interval is the sum of 
HEE of its complementary interval and the horizon length of a static BTZ black hole. This is 


exactly the result suggested in 34 


3.2 Holographic Renyi entropy: Classical part 

In this subsection we develop a systematic way to solve the monodromy problem and calculate 
holographically the classical Renyi entropy for the large interval on a circle at high temperature. 
We need to solve the equation 

= 0, (3.37) 

by tuning the parameters 7 and 6 such that the two solutions for the second order differential 
equation have trivial monodromy along the appointed cycles . Eor convenience, we take a 
conformal transformation 

2ttu 

z = e , (3.38) 
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and define some useful parameters 


ZR = , Z2 = ^ 


(3.39) 


The torus is transformed into a solid annulus with the inside and outside circles being identified. 
The branch cut is at [1, zi\ [j[z 2 , ^], the cycle is the one that goes around the origin for 
n times, and the > 1 cycles are those enclosing [zi,Z 2 ] in different sheets. With the 

conformal transformations, 

z ^ z ■ Zr, (3.40) 

we can cover the full complex plane with a series of annuli. 

To study the monodromy problem for the large interval and high temperature case ana¬ 
lytically, we can take a Laurent or Taylor expansion of the wave function about the origin 
and branch point and also take an expansion with respect to some parameters. Let us first 
consider the expansion about the origin. Because of the monodromy condition, we may rewrite 
the wave function as 

^(i,±) ^ ^|(i±^)j(±)(_2), (3.41) 

where should be single valued in the region Z 2 Zr < \z\ < zi with its Laurent expansion 
being convergent. Assuming the wave function and the parameters can be expanded with 
respect to zr and z as 

OO 

f‘'^\z,ZR,Zi,Z2) = ^fi:^\z,Zi,Z2)zR, (3.42) 

r=0 


OO OO 

7 = ^ JrZR, S = SrZ^R, 
r=0 r=0 

OO 

4^Hz,Zi,Z2) = fr^\zi/^‘2){z/z2r, 

n=—r 


with the normalization as 


(3.43) 

(3.44) 


4 *' = 1 

=0 for r > 0, (3.45) 

we find that and 6r can be solved order by order depending on 7 ^. As zi ~ Z 2 , the 
coefficients /r,n^ are of order 1 and 

z'^rZ^^z^ = (^^)^(^)’-+^, r > 0, n > -r, (3.46) 

Z Z2 

so the wave function should be convergent in the region Z 2 Zr < z < zi. 
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Next we consider the wave function around [ 2 : 1 , ^ 2 ]- The wave function around zi and 
can be written as 


^( 2 ,±) = - Z2)^^^^'^"'^g^^\z,zi,Z2,ZR). (3.47) 

The prefactors encode the information of the monodromy, and is single valued around 
each of the branch points so that it is analytic in the region \z — Z 2 \ < min{ \z 2 ZR — zi\, \ — Z 2 \} ■ 
In this region, the wave function is convergent, and the function can be expanded with 
respect to zr^ {z — Z2) and {zi — Z2), 


00 


g^^'>{z,zi,Z 2 ,ZR) = 

r=0 

(3.48) 

00 

g^^\z,zi,Z 2 ) = ^5(W(2;i,2;2)(2; - Z 2 T 

n=0 

(3.49) 

00 

gi%Hzi,Z2) = gi%\miz2){zi - 2 : 2 )™, 

(3.50) 


m=0 


with normalization 




1 

0 for r > 0 , 


(3.51) 


and the parameters can be expanded as 

00 

70 = X] ")'0,n(^;i - Z2r, 

n=—l 

00 

7r = ^ “ Z 2 T, for r > 0 

n=0 

00 

^ “ Z 2 T■ (3.52) 

n=0 

Taking in the previous result on 5r, we can solve all of the parameters and coefficients order by 
order. We list the solutions of 7 and 6 to the first few lowest orders in Appendix E. Integrating 
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Eqs. (3.13) and (3.17) we get the classical part of nth Renyi entropy 


S„. = 


c ^n+1 c(n+l), , , , ,n 

-ttTR -^--- - log(ni - U 2 ) + c(n + l)(^^(wi - m) 

b n bn obw^ 

TT-^ (n®+ 9n^ - 9n2 + 11 ) 

12960n7 {U 1 -U 2 ) + 

(I9ni° + 85n8 - 125n® + 251n^ - 274n2 + 188) T® . 

- 2449440ni- - 


TT 


TT 


’5878656000ni5 
-126067n2 + 58213)r®(ni - U 2 f) 


(4187n^^ + 12787 n ^2 _ 22521n^° + 55919n® - 100079n® + 143481n'‘ 


4 424 4S44 vr® (I9n4 - 2 n 2 - 11 ) T 6 ^ . 

+c(n-l)(n+l) - 

TfS (4187n® - 572n® - 4218n^ - 1412n2 + 4607) T® g\ 

-12247200»1S- j‘ 

TT® (209n'^ - lOln^ + n^ + ll) 

-^ 20 rf^- 


TT 


r(180041n^° - 104551n® - 17506n® + 43798n'‘ 


'4082400ni3 ' 

-6455n2 - 4607)T®(ni - ^ 2 )^) 6 "^^™ + 0((ni - na)^”) + 0{e 


—67tTR\ 


(3.53) 


Recalling ui — U 2 = I, we find that the classical HRE Sn is in complete agreement with the 


field theory result (2.67) up to order and e 


3.3 Holographic Renyi entropy: 1-loop correction 

In the previous subsection, we calculated the on-shell action for the gravity solution, which 
gives the classical part of the holographic Renyi entropy and entanglement entropy. In this 
section, we derive the 1-loop quantum correction to the holographic Renyi entropy by com¬ 
puting the functional determinants for the fluctuations around the corresponding classical 
background. As proposed in [^, for a handle-body solution realized as quotient space by a 
Schottky group E in pure AdSa gravity, the 1-loop partition function is given by 

00 

Zi-ioop = -Y.Yl - O’ (3.54) 

7 GP 771=2 

±1 

where V denotes the primitive conjugate class of the Schottky group, and q.y ^ denotes the 
eigenvalues of the Schottky group element 7 C SL{2,C), with \qj\ < 1. A group element is 
primitive if it cannot be written as 7 = Z?” for n > 1. 
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To read the 1-loop partition function, we need to find the corresponding SL{2, C) elements 
for every Schottky generator and the primitive elements constructed from them. To study the 
corresponding SL{2,C) elements in the Schottky group, we need to study the monodromy 
around the cycles. We can solve the wave function in different charts covering the Riemann 
surface. If two charts have an overlap, there is an SL{2,C) transformation between the 
solutions in the overlap. For each cycle there are a series of charts covering it so that the 
Schottky group elements are the multiplication of a series of SL(2, C) transformations. The 
crucial point is to study the SL{2, C) transformation between the solutions in two overlapping 
charts. Since we only want to calculate the lowest order terms with respect to the modular 
parameters, we may expand different wave functions in the overlap, such that to each order 
there are only finite number of terms in the expansion. Comparing the coefficients of two 
different expansions, we can read the SL{2, C) transformation to the fixed order. 

Since is expanded at z = 0 and is expanded at z = Z 2 , there is no direct way 

expanding the two wave function in the same region. In the case at hand, to study the Schottky 
transformation we need another wave function, connecting and 0(^0 For convenience, we 
will write the wave function in the new coordinate 

^ Q'n'T{ui+U2) 

^ = (3.55) 

with 

to = —-N-, (3.56) 

g7rT(tt2-'Ul) -|- 1 ' ^ 

which set 

z = 0 —)■ t = —to, 

z = oo —)■ t = to- (3.57) 

The new wave function can be expanded as 

OO - 

0(3)±(0 = 0 + ^ ^ to)z^, (3.58) 

with 

ho{t,to) = '^ho,r{toW, 
r>0 

hm{t,to) = ^ hm,r{to)f , ior m > 0, (3.59) 

r>0,r^2m 

and 

OO 

hm,rito) = ^ ^ hm,r,sto- (3.60) 

s=0 
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Now the convergent region for the expansion is 


l^ol < |i| < |io| 


^—2-kTR _|_ ^'kT{ui—U2) 


^—2'kTR g7rT(ui—^ 2 ) 


(3.61) 


We can solve the wave function order by order with respect to zr and to- 

To study the transformation between and we rewrite in terms of the t 

coordinate 

V,( 2±)(0 = 


(1 + 

•(1 


t-to 


with convergent region 


t |> to- 


q—2itTR _j_ g7rr(«2—111) 


(3.62) 


(3.63) 


^7rT(u2—ui) _ ^—2'kTR' 

By comparing the coefficients in the expansions of two wave functions in the overlapping region 


q—2itTR _j_ g7rr(ii2— 111 ) ^—2-kTR _j_ g7rT(iii— 112 ) 

^0 _ .,,.1 ^1 ^ 1 ^ l^ol 


we get 


where 


g7rT(w2—wi) _ ^—2'kTR 


( 0 ( 2 -) 
0 0 ( 2 +) 


= H 


1 • 


' q—2'kTR _ g7rT('Ui—^ 2 ) 


/ 0(3 ) 


I 0(3+) 


= (-2to)^e5^'^(“i+“2) 


(l+to) 2 ^^“" 01 -to) 2 (l+«^ 

0 


(l+to) 2 (l+s)(l-t(,) 2 (l-s) 


T— 

j,+- 


(3.64) 


(3.65) 


- X-+ 

'j'++ 

(3.66) 


The explicit expression of the matrix elements of T are listed in Appendix F. 

Similarly, we can rewrite the wave function 0(^) in terms of u coordinate and compare 
the expansion coefficients of 0(^) and 0(^) in the region zrZ 2 < \z\ < zi. We read the 
transformation 


where 


/ 0(3 ) 
0 0(3+) 


H 2 — (— 2 fo)' 




2 ■ 



0(3 ) \ 

0(1+) j ■ 

(3.67) 

", ). 

{ZIZ2)~^C^ j 

(3.68) 


The perturbative expansions of the matrix elements C~ and C~^ are listed in Appendix G. 
Because 0(^’+) 0(3’+) (0(i’-) 0(3>“)) share the same monodromy condition around the cycle 
encircling the origin, the transformation matrix H 2 is diagonal. 
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With these wave functions we can get the Schottky generators for the cycle > 1. To 

transform the arguments between different sheets, we need another group element that denotes 
the action of circling around the branch point ui or the origin counterclockwise. Under such 
action, the wave function gets an extra phase so that the transformation matrix is 


M = 


e-2ni 


1 

n 


0 



(3.69) 


With these transformation elements, we can build the Schottky generators for the 
cycle as 

Li = (3.70) 


with i = l,...,n — 1. Ignoring the commutator in the fundamental group, they correspond to 
the thermal cycles in the {i + l)th sheet in the homology group. 

The other Schottky generator corresponds to the horizontal cycle in the first sheet. To find 
the new generator, we need to discuss two other couples of the wave functions to cover the 
cycle. Under a self-mapping conformal transformation 


ZlZ2 

f 

z 


(3.71) 


the energy momentum tensor does not change, so the wave functions under the conformal 
transformation are 


V;(4-±)(5) = V;(l’±)(z)(||)-i = (-l)-i(ziZ2)^^^^^'^»)[l + 


which is convergent in 


I-I A 

Z2 < k < -, 

ZR 


and 




which is convergent at 

Taking a conformal transformation 


z — zi\ < minimi — Z 2 Zr, — — zi}. 

ZR 


(3.72) 


(3.73) 

(3.74) 

(3.75) 


z = 


ZR 


(3.76) 


it is easy to see that the solutions ■ 0 ^^^ and 0 ^^^ share the same convergent region, and 0 ^) 
transforms as 

1 
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Comparing the coefficients of the leading terms, we find 



^ Zj^^'^{-lf^{ziZ2) 2 ^ 

0 ^ 

/ 01 \ 

/ ^( 4 , ) 

V ^ 1 

V 0 

— 1 1 I 

2 ;^" (- 1)2 ( 2 : 12 : 2 )^ / 

Uo J 

^(4,-1-) 


The other wave functions related to by the conformal transformation (3.76) . 

Therefore we get the transformation 


/ ^p^^’ ) \ _ 1 / ) 


Considering the relation (3.74), we have the transformation 

;)(5 

With these results, we obtain the generator 


-) 

+) 


(3.79) 


(3.80) 


Ln = H1H2 




0 


(-1)2 (zi2:2)2- 


0 1 
1 0 


H. 


-1 


0 1 
1 0 


= H, 


{C+)-^C-z^ 

0 


c+(c-)-i4" 


H7^ = HiNHT^ 


h; 


(3.81) 


Up to now, we have built all of the Schottky generators L*, i = 1, • • • re. To calculate the 
1-loop correction to the partition function, we need to find all of the primitive elements up to 
a conjugate. Even though there are infinite primitive conjugate classes, only a finite number 
of them contribute at each order of the expansion with respect to (rei — U 2 ) and In this 
work, we are satisfied to calculate the 1-loop correction up to order (ui — U 2 )‘^ and 

The Schottky group elements can be classified into two classes. In the first class, the group 
elements are generated by 

Li with i < n, (3.82) 


and their inverses. They are similar to the ones in the double interval case as shown in [26| . 
The simplest one is 


7i 


= LiLi+i...Li+k = 


(3.83) 


With this block, all of the group elements in this class are generated by (3.82) and their inverses 
as 

p 

lkuk2,...k,r„m = M^iYl (3.84) 

i=i 

The other class involves L)), which can be written as 


LL = ~ 


1 > 


(3.85) 
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with 


A:' = -log((C+)-'C-z^-) 


(3.86) 


As they take the similar form as (3.84), all the elements take the general form of (3.84). 


However, in (3.84), there are nonprimitive elements, and some of them are conjugate to each 
other. 

In the large interval limit, the asymptotic forms of the group elements are, respectively. 




n n 
'2to 2to 
n n 
2tQ 2to 


., ., / _ H H \ 

iff 1 ~ (^ _|| J 


Then a group element in (3.84 ) has an asymptotic form 

2 ! 

'1k\ ,k2 ,...k2p 


(—rTT( 


h — h 

g n — e ^ " 


) 


—e >1 e n ^ 


— e "n- 


e " 


(3.87) 


(3.88) 


(3.89) 


whose nonzero eigenvalue is 


2 2 p 

n . „ T—r , 27ri 


(^)^n( 

1=2 


±211 h . 

g n _ e „ 


)(1 + O(to))- 


Considering the large interval property for io 


yrT 

^0 ~ ^(^2 - m ), 


we find that the 1 -loop contribution to the partition function is 

OO 

^l-Zoop ^ ^ ^ ^2 ^ ^3 ^ 


(3.90) 


(3.91) 


(3.92) 


s=2 


As the leading order contribution to the 1-loop partition is at order tg^, we only need to 
consider the terms for p < 1, if we are only interested in the result up to order Furthermore 


in (3.84), there may be some Hs that give the same terms as (3.86). Such terms are of order 

_m I , / 

Zpi " , where m is the number in (3.84). Their leading contributions are of order Zr^ . If we 


only consider 1 -loop contributions up to order z^"', we only need to consider m < 1 . 

Here we list the possible primitive conjugate classes, whose contributions to the 1-loop 
partition function are of order no higher than {ui — 1 ( 2 )^ or z^"'. 

1. The group element classes with no L„, including LjLj_|_i...Lj_|_m-i and their inverses. 
Their eigenvalues are 


f 1) - 


2rp2 


vr^T 


1 


= 4n4 sin^ ^ ~ +0{ui-U2y 


(3.93) 
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with degeneracy 2{n — m). Their contributions to the partition function are 


n—1 oo 

Z(i-ioop) ^ ^ 2 (n-m)j;iog(l-((z«r) 

m=l r=2 

= 226800n^ ~ ll)(3n^ + lOn^ + 227)(wi - ^ 2 )^ + 0{ui - U 2 f 

2. The elements and L~^, both of which have eigenvalue 

g(2) = {C+CZ^fz]l'". (3.94) 


The resulting contributions are 


Z 


1—loop _ 


2 ^ 1 og(l-((?( 2 ))^ 


r=2 


r,r 2/nr.. 2( 1 + n^) 2rp2/ \2 

= [1 +- ^^^3 -vr T (ui - U2) 

(-11 - 60n + 15n2 + 12071^ - 9n^ - 60n^ + Sn®) 4 4 


270n7 

+0(2:/?) + 0 ((rti - 712)®)] 

I 3/nr.. ( —l+n^) 2rp2( \2 

+% [1 ^ ^ (“1 “ "“2) 

(-11 - 90n + 15n2 + ISOn^ - 9n^ - 90n^ + hn®) 4 4 


7r"^r^(ui - 772 )"^ 


180n7 


+ 0((77l — 772 )^) + 0 ( 2 ^ 77 )] + 0(zj^)} 

It is clear that the result is expanded with respect to 77i — 772 , zr, and z^'\ 
3. The elements with all kinds of generators include 


7r'‘r'‘(2;i - 22 )^ 

(3.95) 

Ifn 


= (3.96) 


with m = ji — j 2 7 ^ 0. Their eigenvalues are 


<iS> = - + o(.h) + oiT-jT", 

n 

with degeneracy 2n. The contributions to the partition function are, respectively, 


(3.97) 


n—1 00 

Z(l-lor,) _ 2„J/f;iog{l-{,</)’■) 

m=l r=2 

= 2n{^!^(77l- 772)^4'^” —(-251 + 36077 - 11077^+77^) (3.98) 

4'j^4 1 

+ ^^’^^(771 - 772)^44'' —(-1901 + 144077 - 26077^ + 77^), (3.99) 

77* 45 
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for n > 3, and 


for n = 2 , and 


for n = 3. 


r^{l — l00p) 
^2 


4( 


7r"^r^(tti - U 2 Y 

28 


{ZR - 


ry(\-ioop) „(T^^T‘^{nx-U 2 Yi 16 2/3 , 128 

^3 = 6 (-p- 


(3.100) 

(3.101) 


Taking into account all of the contributions, we obtain the 1-loop correction to the nth holo¬ 
graphic Renyi entropy. For n = 2, we have 

- “ 2 )^ “ ^ - “ 2 )^ + ^7r^T^(ui - U2)^^ ZR 

- (2 + ^7r^r^(ni - U2)^ “ "“ 2 )^^ + 0{{ui - U2f) + 0 {zr). 

(3.102) 


For n = 3, we find 


5 . 


1—loop 


16 


177147 


16 


7r"^r"‘(ni - U 2 ) - ' 1 - U 2 ) + 


27 


800 

59049 


TT T (m - n2) z 


2/3 




TT r (m - U2) ) ZR 0{{ui - U2) ) + 0{zY )■ 

(3.103) 


And for n > 3, we obtain 


g\—loop _ _ 


226800n7 
1 


n — 1 


(n -|- l)(n^ -I- ll)(3n^ -|- lOn^ -|- 227)(ni — U 2 Y 

4(n^ — 1 ) 2rr2i 


2 + 


3rfi 


- 7 : T^{ui - U 2 Y 


299 - 444n-M35n224n3 - 3n^ - 12n5n^ 2 /n 

- 2 W-' ^ 

2 (n^ — 1) 2rp2i 


+ 2 + 


n'’ 


-TT T^{ui - U2Y 


6081 - 4626n-h 835n2-h 36n3 - 5n^ - 18n^-h n® 4 ^ 4 , , 4 \ 3 /nl 

-^ ("■ - I 

+0((i.,- «2f) + 0(4''”) (3,104) 


For all the cases, the holographic results are in perfect match with the ones in the field theory 
up to the order we are interested in. 
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4 Conclusion and discussion 


In this work, we completed our study on the Renyi entropy of a large interval on a torus in 
the light of AdS 3 /CFT 2 correspondence. In the case that the interval is not so large, we may 
expand the density matrix in the CFT level by level and compute the entropy perturbatively; 
while on the bulk side, we can follow the prescription in |26] and take into account the size 


dependence 14 to read the holographic Renyi entropy, which is in good agreement with the 
CFT computation. However, when the interval is large, the problem becomes quite difficult. 
On the held side, the perturbative prescription used in the short interval case breaks down, 
and we have to hnd another effective way to compute the partition function. On the bulk side, 
the dual gravitational conhgurations are different from the ones in the short interval case, as 


indicated in the study in 34 


To overcome these difficulties, we developed a new prescription and treatment in both held 
theory and dual gravity. On the held theory side, we proposed in to insert a complete set of 
state bases in the twist sector of orbifold CFT to compute the large interval Renyi entropy. We 
applied this proposal in this paper and focused on the vacuum module of the CFT dual to the 
pure AdSs gravity. We found that the leading linear c contributions were dominated by the 
twist vacuum module and the subleading ones got contributions from all the twist states. This 
allows us to read the leading contributions by applying the Ward identity to the correlation 
function of four twist operators, two at the branch points and the other two at the left and the 
right inhnities of the cylinder. We did hnd the holographic entanglement entropy suggested 
in p4| . 

On the gravity side, we suggested a new set of monodromy conditions on the cycles to 
construct the Schottky generators and corresponding gravitational conhgurations. To check 
the validness of the monodromy condition, we computed the holographic entanglement entropy 
and reproduced successfully the expected value. We read the classical part of the holographic 


Renyi entropy by integrating two differential equations (3.13) and (3.17), one encoding the 
dependence of HRE on the moduli parameter of the Schottky space and the other on the size 
of the torus. Moreover we discussed carefully the 1-loop correction to the HRE, following the 
treatment in [26] . We found good agreements of classical contribution and 1-loop quantum 
correction to the HRE with the leading and subleading large c results in the held theory, up 
to the hrst few orders. For the classical part, the agreement is up to and orders, 

while for the quantum part, the agreement is up to e ^ and orders. 

The study in this work presents another piece of evidence to strongly support the holo¬ 
graphic computation of the entanglement entropy in the context of the AdS 3 /CFT 2 corre¬ 
spondence. Taking into account the accumulated evidence on the holographic Renyi entropy 
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in the cases including double-interval and single short interval on the torus, it suggests that 
the holographic computation is exact perturbatively not only at the classical level but also at 
thel-loop quantum level. Furthermore, our field theory study shows that there are actually 
1/c corrections in the partition function when the Riemann surface is of higher genus than 1. 


It would be interesting to see if the agreement could go beyond the 1-loop level 28, ^ or even 
nonperturbatively. 

Our study could be generalized to other cases. In particular, it is interesting to study the 
higher spin Renyi entropy of the single interval on a torus by direct field theory computation 


42l-l46|and Wilson line prescription in the bulk 147-49 


The study of holographic entanglement entropy may shed light on the AdSs quantum 


gravity 24,50 . There are two essential questions on the quantum AdSa gravity. One is on 
the precise definition of the quantum gravity, string theory, or something else. The other is 
on the construction of the dual CFT. There is ample evidence, for example, the work in |24| , 
that the dual CFT might not exist. However, the results in this work and other related ones 
suggest that there exists an equivalence between the semiclassical AdSs gravity including the 
pure gravity sector and the large central charge limit of a 2D CFT, which has a sparse light 


spectrum 21,23 51-^. In our study, it turns out that in the large central charge limit, the 
vacuum module dominates the contribution to the partition function. It is not clear when 
the states in the sparse light spectrum begin to contribute. Moreover, the regulated on-shell 
action of the gravitational configuration in the AdSs gravity is a Liouville theory. This raises 
the issue if the dual CFT could be a Liouville CFT. For a recent study on this issue, see [54| . 
It would be interesting to see if it is possible to prove the equivalence by using the Liouville 
theory]^ 
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Appendix A States in the vacuum module 

In this section, we list some low-lying excited states in the vacuum module. We focus only on 
the holomorphic sector. For the antiholomorphic sector, it is similar to the holomorphic one. 
^We would like to thank the anonymous referee for pointing out this possibility. 
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The first few excited states up to level 4 are, respectively, 


2 ) = y-T_ 2 | 0 ) 


3) = 


I 4,2) = (| + He)-|(L_2L_2-^L_4)|0). (A.l) 

The corresponding vertex operators at the origin and the infinity take the forms, respectively, 

L -2 T{w) |^„=o 

—5- w^T{w) I 

W^OQ 

L _3 dT{w) 1^=0 

-w^dT{w) - Aw^T{w) l^^oo 
L _4 —)• -d‘^T(w) |to=o 

^w^d‘^T{w) + 5w'^dT{w) + 10w^T{w) |^^.oo 

L_2T_2-^L-4 ^ :T{wf :-^d^T{w)\^=o 

5 lU 

w^{: T{wf : -^5^r(u;)) \w^oo 

Appendix B Conformal transformation for : T‘^{z) : 


In the calculation, we need the conformal transformation of : T‘^{z) :, which is not a primary 
operator. Under a conformal transformation z —)• w{z), we have 

N X ,dW ,0 C r ^ . X 


where 


Tiz) = nw)i'^f + ^{w,z}, 


u;"'3 tc" 


is the Schwarzian derivative and the prime denotes the derivative with respect to z. For 
: T‘^{z) :, we have 


■■T\z)-. = /^^T(zi)T(z) 

J zi- z 


,wr^w"' 


w" w"" 


——{-22hi — T + 480( —)^— - 100( —- 180— — 
1440 ^ ^ w' ^ w' w' ^ w' w' w 

1 in" 2 in'" 3 in" 

+36^) + 2T{w){w'f{-\C^f + \^) + dT{w){w'fl^, 

w 4: w 6 w 2 w 

+ : T\w) : {wr + rM(u;')^ 4 + 4)^- 
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Especially for the conformal transformation 


— V 

w = , 


(B.4) 


we have 




C, 


(B.5) 


11 c, Ti^ , c 10 ,,,( 9 ui ,9 tt^ ~.,dw-,n vr ,9 




'36 ' 90 ^n4/34 ' ^ 3 ' 3 

Appendix C Correlation functions /, 




In this computation of the Renyi entropy by inserting the twist sector states, we need to 
compute the correlation functions /*. Here we list the results for the first few ones needed in 


the relation (2.32) 


log(f I E(0) I f) = c(- 


(—1 + n^) 7 r^ 2 (“11 + 20n^ — 18 n^ + 8n^ + n®) 7 r^ 4 


4^ + 


+ 


72,d%3 ” ' 25920/34n7 

(-2497 + 2160n2 + 294n^ + 40n® + 3 n®) 7 r^ 4 


n 


453600/3%7 


r + o(p), (c.i) 


/2 = (1 + 


2(-l + n 2 ) 7 r 2 2 (299 - 444n + 135n2 + 24n3 - 3n^ - 12n'^ + n®) 7 r^ 4 


3/3^n^ 


- 


1 (733 - 672n2 - 63n^ + 2n^)'K‘^ . 5 


= 1 for n = 1 


945/3%7 

^2 


54/3^n^ 

+ O(l^) for n / 1, 2 


r) 


TT" ,9 257r^ 


TT^ 

= c—_+ 11 + +_ 

128,04 ‘ 1 ^4/32 ^ 768,04 


1 37r4 

+-+ 0{l^) for n = 2, 

^ c 128,04 ^ ^ 


(C.2) 


/3 = (1 + 


(-1 + n2)7r2 2 (6081 - 4626n + 835n2 + 36n3 - 5^4 - 18^^ + n®)7r4 . 


-P - 


/3‘^n^ 

1 (-933 + 1148n2 - 217n4 + 2n®)7r4 


36l3^n'^ 

+ 0{l^) for n / 1,3 


630/3^n'^ 


n 


1 for n = 1 
1287r4 


6561,0^ 


^ (1 0~ ++09^ 


87r2 ,2 2807r4 , 1 327r4 


27/32 6561/3' 


-D + - 


c 6561/34 


/^ + 0(/^) forn = 3, (C.3) 
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/4,1 = (1 + 


4(-l + n^)'iT‘^ 2 (201919 - 108120n + 14115n2 + 24071 ^ - 39n^ - I20n^ + 5 n®) 7 r^ . 


3f3^n^ 


-V - 


1 2(83213 - 136080n + 55020n2 - 2163n^ + 10 n 6 ) 7 r 4 . . 


4725/34n7 


135/34n7 

Z^ + 0(Z'’) forn^l,2,4 


1 for n = 1 
Svr^ 

'64/^4 


,4 TT^ ,2 4:737r‘ 

c^Z" + (1 + 2;^Z2 + 


1920/3^ 


1 37r^ 
c 320/34 ‘ 


^ +-^, 7 , 7754 ^^+ 0 (Z'^) forn = 2 


1257r4 


Svr^ 


-I + (1 + 2 ^ 4“ 


226l7r4 


'4096/34 ^ 16/32 49152/34 


n + 7 


1 997r^ 


c 40960/3' 


-l^ + 0{l^) for n = 4, 




(C.4) 


/4,2 — (1 + 


4(-l + n 2 ) 7 r 2 2 (299 - 456n + 135n2 + 48n3 - 3^4 - 24n® + n^)-K'^ . 




-V - 


27/34n7 


n 


1 2(124247 - 166320n + 43260n2 - 1197n^ + 10 n 6 ) 7 r 4 . . 


1 for n = 1 


4725^4^7 

3377r^ 


r + 0{r) forn/ 1,2 


+ + + _ 

64/34 ^2/32 ^192O;04 


+ forn = 2 . 

^ c 32 O ;04 


(C.5) 


Appendix D Mutual Renyi information for the double inter¬ 
vals 


In this appendix, we list the Renyi mutual information for the double intervals, which has 
been computed in [^. In terms of a small cross ratio x, the leading and next-to-leading 
contributions are, respectively, . 


rtree 


{x,x) = 

c{n — l){n + l)^x^ ^ c(n — l)(n -|- l)^x^ ^ c{n — l)(n -|- l)2(1309n^ — 2 n^ — ll)x^ 


288n3 288n3 414720n7 

c(n - l)(n -F l)2(589n4 - 2n‘^ - ll)x^ 

207360n7 

c(n - l)(n -F l)2(805139re^ - 4244n® - 23397n^ - 86n2 -h 188)x® 

3135283207^1 


-|- 0(x’^) -|- c.c(D.l) 


{n + l){n‘^+ ll){3n‘^+ Wn'^+ 227)x^ 

7257600n7 

^ (n -F 1)(109778 -F 1495n® -h 11307n4 81905772 - 8416)x5 

119750400n9 

, (77 -F 1) (1444050771° -h 1911297477® -h 14056530577° -F 100052783777^ - 167731255772 - 14142911)x° 

10461394944007711 

-|- 0 (x’^)-|-c.c. (D. 2 ) 
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Appendix E The accessory parameters 


After imposing the monodromy condition, the accessory parameters 7 and 5 can be solved 
order by order. Here we just list the expansion coefficients of the first few orders 


70,-1 = 


27r(n — l)(n + l)Tz 2 


\ 2 f^ I 


■) 7o,o = 


7 r(n — l)(n + 1)T 


7o,i = 


7 r(n — l)^(n + lyT 


I ^\2„ 


U 70,2 = 


(—1 + n)^(l + nfirT 


70,3 = 


70,4 = 


6n^Z2 

7 r(n - l)2(n + l)^ (229n^ - 2n^ - ll) T 


4320n8z| 


7 r(n — l)^(n + 1)^ (l09n^ — 2n^ — ll) T 


2880n8z| 


70,5 = 


70,6 = 


7 r(n - l)2(n + 1 )^ (6299911® - 1724n® - 953711^ - 86n‘^ + 188) T 


21772801112^5 


7 r(n - l)2(ii + 1)2 (2015911® - 716ii® - 3993ii^ - 86n^ + 188) T 


870912ni2zf 


7r(n — l)2(ii + 1)2 t 
“ 15676416000ii164 


(299803787: 


n 


12 


- 1284030611 ^° - 7225394711® - 2555548ii° + 5564373n'^ + 67854ii2 - 58213) 


71,0 = 0, 71,1 = 0, 71,2 = 0, 

7 r(ii — l)®(n + 1)®T 7 r(ii — l)®(ii + 1)®T 

=-SAi-• =--’ 

7 r(n — l)®(ii + 1)® (439n^ — 2ii^ — ll) T 

“ 432 O 11 I 04 ’ 

7 r(ii + 1)® (I99nl’ - 597ii° + 595ii® - 193n^ - 17n® + 35n2 - 33n + ll) T 
~ 1728ni04 ’ 

■JT ( Yi, I 1 ^ ^^7”' 

- ^ 4, - (4110527n®^ - 12331581n®° + 12268009n° - 3919811n® 

32659200ni44 ^ 

-541434n^ + 1115726ii° - 1053566ii® + 354954n'^ + 371n® - 12409n2 + 13821ii - 4607) 


72,0 = 0, 

7r(ii — 

72,3 = - 


72,1 = 0 , 72,2 = 0 , 

1)2(11+1)2 (3ll2-l)r 

6n®4 
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72,7 


72,4 = 


7r(n — l)^(n + 1)^ (3n^ — l) T 


Arfiz2 


72,5 = 


72,6 = 


7r(n — l)^(n + 1)^ (l469n® — 521n"^ + + ll) T 


1440nio4 


7r(n — l)^(n + 1)^ (749n® — 281n^ + v? + ll) T 


576nio,z6 


7r(n + 1) T /';^7287061n^2 - 34574122n^i + 10493170n^° + 13587782n® - 6779897n® 

10886400ni44 V 

-27988n^ + 404292n® - 780596n^ + 383843n'‘ + 12910n^ - 1106277^ + 9214n - 4607) 


2Tr^T‘^ 

<5o =-2—’ = <^1,1 = <^1,2 = <^1,3 = 0, 




7r^(n^-l)^T2 7r^(n - l)^(n + l)^r^ 

-36rfVl-' =- 


18n®^ 


TT^ (n2 - 1)^ (I039n^ - 2n2 - ll) 


<Jl,6 = 


12960nio4 


<^1,7 — — 


77^ (n2 - 1)^ (43977^ - 277^ - 11) r2 


432O77IO4 


<^1,8 = 


7r2(n - 1)3(77 + 1)3 (1572646777® - 11649277® - 64177877^ - 1412772 + 4607) 


1306368007714^1 


(^1,9 = — 


77^(77 - 1)3(77 + 1)3 (447214777® - 4593277® - 25369877^ - 1412772 + 4607) 

326592OO77I44 


61 10 = w ( 713275944777^2 _ 9310644677^® 

47029248000771®^® ^ 

-51582518777® - 492990877® + 1605783377^ + 73554772 - 78893) 


<^2,0 — <^2,1 — <^2,2 — <^2,3 — 0, 

^2 i^2 fQ^2 ^\rr2 

(^2,4 = 


TT" (772 — 1) (3772 — 1) T2 


677®4 

^2/^2 /o^2 


^2,5 = - 


772 (772 - 1)" ( 3 n 2 - 1) T2 


377® Zn 


7r2 (772 - 1)^ (326977® - II2I774 + 772 + 11) r2 


<^2,6 = 


2I6O77IO4 


(52,7 = — 


7r2 (772 — 1) (146977® — 52I774 + 772 + 11) T2 

72077104 
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TT^ {p? - l)^ (56156801ni° - 20579551n® + 40454n® + 681358n^ - 6455n2 - 4607) 

“ 21772800ni4zf 

7r2 (n^ - l)^ (I7031281n^° - 6432271n® + 5174n6 + 293278n^ - 6455n2 - 4607) 

“ 5443200^14^9 

7^2 (n‘^ - if t2 

(52 10 = - 1(28908872157n^^ - 11219722453n^2 

7838208000ni8zio ^ 

-16128831ni° + 654870399n® - 21992209n® - 17014839n^ + 18888371^ + 78893) 


Collecting all these coefficients and changing back to the u coordinate, we find 


70 = (4-1) 


1 


+ (n^ - 1) TT^r^ |(-^)(ni - U2) 


(n® -|- 971.4 _ _|_ 215 

+(-giSis-><“' - 

, 7r4r4 (I9ni° -b 85n® - 125ni^ -b 251n4 - 274n2 -b 188) _ . 

+ (-^1 - 112 )' 


68040ni2 


+ 


776706 


(4187n^^ -b 12787n^2 - 22521n^° -b 55919n® - 100079n® 


122472000ni6 

-bl43481n^ - 12606777^ -b 58213) (tti - 772)^ + 0((7ti - 172 )®)} 


71 = 7 rT (77 - 1 ) < 7777 ( 771 - 772 ) -b 


TT^r^ (19774 - 277^ - 11) 


9776! 1357710 

774704 (41877^8 _ 572^6 _ 42i8n4 - 1412772 -b 4607) 
25515 O 77 I 4 


(m - 772)^-b 


(771 - U2f + 0{{UI - U2f) 


72 = TT T^{n - 1) 


2 I 4 (3772 - 1) 


377® 


-(771 - U 2 Y 


7r2r2 (20977® -10l774 -b 772 -b 11), 

+---1(“0 - 

7r4r4 (I8OO4I77IO - 10455177® - 1750677® -b 43798774 - 6455772 - 4607) , 

+-^. .. .^(771 - 772 )^ 


-bO((77l - 772 )®)} 


850507714 


<5i = 


- 0’ If b~ (“1 - 


+ 


4057710 

c 

7r4T4 (418777® - 57277® - 4218774 - 1412772 -b 4607) 


-bO((77l - 772 )®)} 


5103007714 


(771 - U2f 


49 



= 7r^r®(n2-l) 


2/8 (3n2 - 1) 


-{ui - U 2 ) 


47r2T2 (209n® - lOln^ + + ll) ^ . 

+-^ 

vr^T^ (I80041ni° - 104551n® - ITSOGn^^ + 43798n^ - 6455^^ - 4607) , 

+-^-8505^^-i(“. - “2)' 

+0{{ui - lf2)®)} 

Appendix F T matrix 

In this section, we list the leading order terms in the expansion of the T matrix in the Schottky 
transformation. For any matrix element, we may expand it as 

00 00 

T = Y.Y1 (F.l) 

r=0 n=—1 

There are relations among the matrix elements 

T-rt = i-^TT-- (F.2) 

T++=T—\n^-n (F.3) 

r+-= r-+L,^_n fF.4) 


For the matrix element T , its expansion coefficients are, respectively. 

To”! = ^0,0=-^ 


36^3 

_ —+ 6n^ + — 6n + 4 

36n® 

-547n® + 3300n'^ + 1450n6 - 432071® + 489n^ + 900n® - 2600n2 + 120n + 1208 

32400n7 

2563n® - 4500n^ + 230n® + 4320n® - 452177*^ + 270077® - 28077^ - 252077 + 2008 

32400779 


l * — 677® — 577 ^ + 677 + 4 


^0,2 — 


^0,3 — 


^0,4 ~ 


285768007711 


(35130177®^ - 215413877^^ - 48272777®° + 302143877° - 132016577® - 44381477®” 


+216517977° - 76704677® + 20269277® + 32575277® - 151855277° + 1780877 + 602272) 


Tq q = 157309777®° + 333484277®® - 63640577®° - 302614277° + 400841777® - 325823477®” 

+92966577° + 199949477® - 285322877® + 243079277® - 110796077° - 148075277 + 1232608) 


50 



^ 0,7 


-^(-50020417977^® + 3117021000n^® + 239603780n^^ - 447758664077^3 

51438240000773 ®^ 

+289890280677^2 + 23660280077^^ - 3339143740773° + 192280992077° - 160924861977® 
-247284600773' + 274449064077® - 88942392077® + 122304353677^ + 32578080077® - 2411735680772 
+1208064077 + 754291456) 


^ 1,-1 ~ ^ 1,0 ~ ^ 1,2 ~ 0 ) 


^iT = 


4 (772 — 1) 
377 


ny = (4(12 + 577 - 23772 - 1077® - 277^^ + 577® + 1377®))/(4577®) 


T’m = - 


2 (1377® - 1277® - 2477^ + 2477® + 9772 - 1277 + 2) 


27it7 


rrr =-—^(-101693773® - 11676077° - 3331777® + 27552O773' + 26736177® - 20160077® 

3’® 8505077° 

-5883177^ + 4368077® - I3OI68772 - 84077 + 56648) 

_ __ _ _ 4 (772 — 1) 

^2,-1 = ^2,0 ~ ^2,2 = *3) ^2,1 = 

__ _ 4 (6277® + 7577® + 1477^ - 9077® - 112772 + 1577 + 36) 

2,3 4577 ® 


Appendix G C matrix 


For the matrix elements C~^ and C , their expansions are similar 

00 00 

e- = EE+«*«*S (G'l) 

r=0 n=0 


and 


C+ 


C- 


n^—n 


(G.2) 


Here we list the ones for C : 


^ 0,1 ~ ^ 0,3 ~ ^ 0,5 ~ ^ 0,7 ~ 0 


ri- ri- — ~ + 77 + 2 

^ 0,0 - ^ 0,2 - ya 

— 67773' — 14077® + 11477® + 48773 — 6377® + I8O772 + I677 — 88 

108077^ 


C, 


0,6 


-T-(-49197733 - 10570773° + 1068177° - 32277® - 6053773’ 17540^6 

1360807731 ^ 

-289377® + 2506773 + 387277® - 15176772 - 68877 + 6016) 


51 



c, 


0,8 


-^(-6094181n^^ - 13347880^^"^ + 15437132n^^ - 4050896n^2 - 8054742n^^ 

24494400071^5 v 

+2450184077^° - 860627677° + 1316172877® + 513262777'^ - 2338724077° + 590505677° - 1337956877^^ 
-441670477° + 2395328077^ + 69708877 - 7451264) 


^ 1,0 ~ ^ 1,1 ~ ^ 1,3 ~ ^ 1,5 ~ 


4 ( 77 ^ — 1 ) 

377^ 


2 (8377° - 9077° - 4277^ + 18077° - 93772 - 9077 + 52) 

135776 


c 


1,6 


—-^(3637777^° - 4225277° - 3653577® + 7249277^ - 1868577° + 277277° + 1225977"^ 

1701077^0 

-5401277° + 1328877^ + 2100077 - 6704) 


C2A = 


_ _ _ _ 4 (772 — 1) 

^2,0 ~ ^2,1 = ^2,3 ~ ^2,2 = ^^2 

2 (14977° - 9077° - 10277^ + 18077° - 99772 - 9077 + 52) 

45776 
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